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Abstract 

Infrared observations at an effective wavelength 
of 3.5 vi of a selected group of Long Period Variable 
(LPV) stars are presented. Mira type and semi- 
regular stars of M, S, and C spectral classifications 
are monitored throughout the full cycle of variability. 
Although the variable infrared radiation does not 
exactly repeat in intensity or time, the regularity 
is sufficient to produce mean 3.5y light curves. 

For the M, S, and C stars as a group, the mean 3.5y 
alterations are about 0^7 . The 3 . 5y maximum 
radiation lags the visual maximum by about one- 
seventh of a cycle, while the minimum 3.5y intensity 
occurs nearly one-half cycle after infrared maximum. 

In some stars, there are inflections or humps on the 
ascending portion of the 3.5y light curve which may 
also be seen in the visual variations. 

The 1 . 04y narrow band photometry of Lockwood 
and Wing (19 71) and Lockwood (1972) , in conjunction 
with the present 3.5y observations, is utilized to 
form a color index ( [1. 04y] - [3 . 5y] ) . This color 
index is available during most of the full variability 
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cycle for the stars in common between the two programs . 
The color index is converted into a color temperature 
under the assumption that a Lyra (AOV) is zero magnitude 
at all wavelengths and, in the infrared, radiates as 
a black body with a 10,000 °K color temperature. 

From the deduced color temperature, via the 
Planck function and the 3.5y magnitudes converted 
into absolute fluxes, the angular size of the source 
is calculated. Color temperatures, sizes, luminos- 
ities, visual magnitudes and radial velocities are 
computed for the cool M giant variables over the 
entire cycle and compared with observational results . 

The color temperatures attain maximum coin- 
cidentially with visual maximum or zero phase, while 
the minimum is realized almost exactly one-half 
cycle later. Maximum compression occurs about one- tenth 
cycle before maximum temperature; maximum radius 
occurs about four-tenths of a cycle after minimum 
radius. Maximum luminosity happens about one-tenth 
cycle after maximum temperature and minimum energy 
output takes place nearly one-half cycle after 
maximum. The fractional change of temperature is 
about 20%; that of diameter is a little larger and 
the luminosity alters by a factor of two. 

The calculated results are in excellent agreement 
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with Pettit and Nicholson (1933) and in good agreement 
with other observed quantities. The phase lag of 
maximum luminosity to minimum radius is similar to 
that seen in Cepheid variables and suggests that the 
physical mechanisms producing the variations may 
be similar. 


Tr 
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...“Variability is a disease of bloated stars"... 

P. W. Merrill, 1938 
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I 

INTRODUCTION 

One of the most curious facts of history is that 
the first variable star ever discovered is still 
probably one of the least well understood. A second 
point is that, even with the great interest in the 
heavens demonstrated by early civilizations, the first 
recording of a star altering its apparent brightness , 
other than supernovae, waited until 1596. On August 13, 
1596, a German clergyman, David Fabricus , announced that 
the star designated omicron in the constellation of 
Cetus , the whale, did not appear constant in intensity. 
Again in 1638, yet fourteen years before the invention 
of the telescope, this same star, omicron Ceti , was 
rediscovered by Holwarda and its variability was 
confirmed. This interesting star, o Ceti, was given 
the name Mira - the wonderful - in honor of its unusual 
behavior. 

Several theories were formulated to explain the 
underlying cause or causes of the stellar variability. 
One of the earliest and most obvious hypotheses was a 
binary or multiple star system with eclipses producing 
the changes in visual magnitude. This theory is able 
to account for many of the variable stars -- namely 
the extrinsic variables. 
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Some brightness alterations, however, were not 
acceptably explained by the binary star hypothesis. 
Shapely pointed out that the probable dimensions of 
some of the variable stars were on the order of ten 
times larger than the computed orbits of their 
companions. Gradually, as theoretical work on stellar 
structure progressed, the pulsation hypothesis, as 
supported by Eddington, was invoked to account for the 
stellar variability. At the present time, it is 
generally accepted that many stellar luminosity changes 
may be adequately described by radial oscillations 
or pulsations of either the entire star or its 
atmosphere alone. These stars are then called 
intrinsic variables. 

One particular type of variable star called 
Cepheids , after the prototype of this class 6 Cephei , 
has been quite extensively studied. From an observa- 
tional standpoint, these stars are relatively bright, 
have well defined periods, and are quite reproducible 
from cycle to cycle. Due to the accuracy and availa- 
bility of the observational data, many theoretical 
models have been proposed. Recent work has quite 
satisfactorily explained the static structure of F, G, 
and K giants from which dynamical calculations are able 
to proceed. At the present time, the observable 
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parameters of Cepheids (apparent brightness, radial 
velocity, and period) and the calculated quantities 
of Cepheids (temperature, angular diameter, and 
luminosity) are quite satisfactorily explained by 
the stellar structure dynamical models. 

On the negative side, however, there exists a 
group of intrinsic variables which are not as well 
understood as the Cepheids . These are the red giant 
stars called Mira variables after the first one of this 
type discovered. Mira type stars are part of a larger 
group designated long period variable (LPV) stars. 

The Miras are characterized by cool photospheric temper- 
atures (~2500°K) , large size (~100 R Q ) , and periods 
approximating one year. 

Long period variables have been studied since the 
turn of the century and have been quite extensively 
described and catalogued. However, the physical 
processes which cause and maintain the pulsations are 
less assured. Due to the cool effective temperature 
of the photosphere, the maximum energy is radiated at 
infrared wavelengths longer than about one micron. 

Also, due to the low temperatures, LPV stars possess 
extremely complex spectra. Absorption by atomic lines 
and molecular bands and line emission by hydrogen distort 
the energy distribution. As a result of these complex 



4 


spectra and convective energy transport, theoretical 
models of non-variable red giants have only recently 
been forthcoming, whereas dynamical models for the 
variability have been extremely scarce. Because of 
this dearth of prediction, comparison with observa- 
tional quantities has been non-existent. 

The original goal of this research program was to 
ascertain if stars do alter in intensity in the infrared 
at a wavelength of 3 . 5y . If so, the aim was to deter- 
mine the magnitude of the variation and its phase with 
respect to the well documented visual brightness 
changes. The wavelength of 3 . 5u is relatively unaffected 
by interstellar and circumstellar extinction, at least 
compared to the visual region, and cool stars radiate 
much of their energy at this wavelength. Because the 
period and amplitude of Mira type stars are more 
reproducible than in other types of cool giants, the 
decision was made to monitor a selected group of Miras 
to ascertain answers to the previously posed problems. 
Also included in the observational program are a few 
semi-regular variables to determine some of the same 
quantities . 

The general tone of this present investigation 

be primarily the presentation of observational 
results. Quantities such as color temperature, angular 
size, luminosity, radial velocity and others, which 
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may be calculated from the observational data, are 
included in later chapters. It will not be the attempt 
of this essay to explain all of the phenomena which 
occur in the colors and spectra of variable red giants. 
It is hoped, however, that the quantities and relations 
presented herein will inspire further development of 
stellar models to explain the complex nature of long 
period variable stars. 
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II 

CHARACTERISTICS OF LONG PERIOD VARIABLES 

The ensuing paragraphs summarize a few of the 
more salient features of long period variable stars 
and are included for historical, descriptive, and 
informative purposes. 

Periods : The frequency of occurrence of LPV's 

maximizes near periods of 280 days (Merrill, 
1938) with the later types possessing longer 
periods (Cameron and Nassau, 1956) . One of 
the outstanding differences between LPV stars 
and the Cepheids is that the cooler giants are 
not nearly as regular in period (-10%) as the 
shorter period, hotter stars. These single 
cycle alterations of period appear random and 
there is no systematic change. There exist two 
exceptions to this rule, namely R Aql and R Hya, 
both whose periods have decreased about 20% 
(Campbell and Jacchia, 1941). A period-spectrum 
relation does exist for LPV's with the later 
spectral types having longer periods (Campbell 
and Jacchia, 1941) , but there seems to be no 
period-luminosity law (Glasby, 1969) . 

Spectral Type : Approximately 89% of the long 

period variables are classified as M stars, 7% 
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as R, N and C stars and the remaining 4% as 
S type stars (Merrill, 1960) . The oxygen rich 
stars are denoted as M; the carbon rich are 
given the R, N or C classification, whereas the 
S type stars exhibit relatively high heavy 
metal abundances (Merrill, 1960) . Eggen (1972) 
states that all N type stars are variable and R 
types are constant with a few exceptions . 

Wyckoff (1970) finds that the average change in 
spectral class from maximum to minimum is three 
subclasses for the M type stars. 

Galactic Distribution and Age : In general , long 

period variables are thought to be old, evolved 
objects. None are found in Pop I galactic 
clusters but few are found in globular clusters 
(Payne-Gaposchkin, 1954). The shorter period 
Miras have higher space velocities (Eggen, 1955) 
and tend to form a spherical distribution 
throughout our galaxy, whereas those of longer 
period tend towards an intermediate population 
(Payne-Gaposchkin, 1954) with lower space 
velocities and somewhat higher masses (Glasby , 
1969). In contrast, the non-variable M giants, 
as well as the R and N stars, tend toward a 
rather flattened galactic distribution (Payne- 
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Gaposchkin, 19 54) . The latest M type Miras of 
the longest periods may even be Pop I objects 
(Cameron and Nassau, 1956; Eggen, 1955). 

Visual Light Curves : The amount of the apparent 

visual brightness alterations (Am v ) increases as 
the period lengthens. The mean change in m^ is 
4 m for variables with periods near 200 days and 
Am^ increases to 5 m near 300 day periods 
(Merrill, 1960b) but tends to remain the same 
for longer periods (Merrill, 1960a). Mira type 
variables often have a pronounced hump either on 
the ascending or descending portion of the light 
curve which may not be present during each cycle 
(Glasby , 1969). S type variables generally 
possess wider flatter maxima and a greater visual 
magnitude alteration than do the M type stars 
(Campbell and Jacchia, 1941) . R stars have broad 
flat maxima with sharp narrow minima; N type stars 
evidence smaller Am^, longer periods, broad maxima 
and narrow minima, decline faster than ascend 
and possess humps as a rule rather than as the 
exception (Campbell and Jacchia, 1941). Terrill 
(1969) demonstrated that the M type Mira variables 
show a tendency to be somewhat brighter at a given 
spectral type on the rising portion of the light 
curve than on the declining branch. 
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Spectra ; The variables designated as Mira type 
are characterized by hydrogen emission lines, 
whereas cool variables , in which the lines are 
absent, are called semi-regular (SR) or irregular 
(Irr) (Campbell and Jacchia, 1941) . These bright 
emission lines appear approximately one-quarter 
cycle before visual maximum and persist for about 
one-half cycle until about phase 0.25 and then are 
absent until the next cycle (Campbell and Jacchia, 
1941) . These hydrogen lines are believed to 
arise from shock waves in the stellar atmosphere. 
The shock wave is generated below the reversing 
layer within the hydrogen convection zone and 
then propagates outward through the reversing 
layer (Glasby, 1969) . 

Luminosities : Clayton and Feast (1969) have 

determined the absolute visual magnitude at mean 
maximum light for the M type Miras. They 
reported an absolute magnitude of -3 m for Miras 
of 200 day periods and decreasing to -l m for 
the 400 day periods. They also report the 
absolute visual magnitude at mean light is 
about 2™ dimmer than at mean maximum. The Mira 
type S stars typically have M about ~3 m , also 
(Takayanagi , 1960). The absolute bolometric 
magnitude for LPV stars as a group is about -4 m 
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(Campbell and Jacchia, 1941) but near -5 m for 
the Miras only (Smak, 1966). In M type stars, 
the change in visual blanketing throughout the 
cycle, caused by molecular bands attributed to 
TiO, could account for some of the differences 
as compared to stars of other spectral types and 
also could explain the differences between the 
visual and bolometric light curves (Campbell 
and Jacchia, 1941) . 

Masses : Mira type variables were believed to 

have masses on the order of ten solar masses 
but recent work by Fernie (1959) and Fernie 
and Brooker (1961) on X Oph and o Cet seems 
to indicate masses more on the order of one 
solar mass. 

Mass Loss : Many of the brighter M type Mira 

variables exhibit excess radiation near lly 
(Gillett, Merrill and Stein, 1971; Gehrz and 
Woolf, 1971) which has been attributed to 
circumstellar shells composed of silicate grains 
(Woolf and Ney, 1969). Gehrz and Woolf (1971) 
deduce an average mass loss for Miras of about 
2 x 10 ^ M^/year. 

At the conclusion of this recap of characteristics 
of long period variable stars , a short summary of the 
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very early yet often cited work of Pettit and 
Nicholson (1933) on LPV stars is included. They 
measured eleven stars throughout the variability cycle 
and found an average radiometric magnitude change of 
nearly one magnitude with the maximum energy output 
occurring about one-seventh cycle after visual 
maximum. They also computed the stellar diameter to 
be more in phase with the luminosity or about 180° 
out of phase with the temperature. Their deduced 
mean diameters were 18% larger at minimum than at 
maximum light. 
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III 

EQUIPMENT 

All of the previously unpublished observations 
were obtained with detectors , dewars , and photometers 
designed and built at the University of Minnesota 
under the direction of Dr. E. P. Ney. 

Nearly all of the observational data which forms 
the basis of this thesis was obtained with the 30 inch 
(76 cm) Cassegrain telescope of the O'Brien Observatory 
of the University of Minnesota located near Marine-on- 
St. Croix, Minnesota. A few data points were 
accomplished via the 60 inch (152 cm) metal mirror 
telescope of the University of Minnesota - University 
of California, San Diego Infrared Observatory on 
Mt. Lemmon near Tucson, Arizona. 

In general, throughout the time span of June, 

1968 through June, 1972, when the data were obtained, 
all equipment was altered and, hopefully, improved. 
Assortments of detectors, dewars, beam switching 
techniques, electronics, interference filters, and 
observers were utilized. One of the prime concerns 
during the system alterations was to maintain, as 
nearly as possible, the same effective 3.5y wavelength 
of all photometric systems, so that meaningful inter- 
comparisons among all the data could be accomplished. 



Figure 1 is a block diagram of the early (1968) 
single wavelength photometric system. 

From 1968 through June, 1971, the beam switching 
system consisted of a four section, two level mirror 
which, when rotated at 5 Hz, produced a 10 Hz signal 
on the detector. Two sections of the mirror were 
movable, allowing the separation of beam centers to 
be varied from 1 mm. through 1 cm. corresponding to 
25 Ti through 25o7l , respectively, with the scale of 
the O’Brien telescope. 

Most of the observations made in Minnesota 
utilized a detector system employing a 3 mm. focal 
plane limiting aperture which resulted in a 1-1/4 
arc min beam diameter. The rotating two level mirror 
was generally used with a 2 mm separation implying 
that the two beams overlapped in area on the sky. 

This overlap introduced no problems on stellar point 
sources and needed no correction factors. 

In June of 1972, the rotating mirror was replaced 
by a linearly translating mirror. The cam driving 
mechanism produced a 10 Hz square wave signal on the 
detector. This translator had a fixed beam separation 
of 1-1/4 mm and only systems with focal plane dia- 
phragms less than or equal to 1 mm were used. 

The first infrared (3.5u) observations obtained 
with the Minnesota equipment were made in January, 



TELESCOPE 



Figure 1. Block diagram of the equipment. 
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1968. These earliest detectors possessed a Noise 

-13 

Equivalent Power (NEP) of about 10 watts in a 1 Hz 
bandpass. This sensitivity and resulting signal-to- 
noise ratio imposed some restrictions on the possible 
choices of program stars, more of which will be mentioned 
later. As systems improved after 1968, the NEP was re- 
duced by about a factor of 10. Thus, the accuracy and 
rate of observations improved, allowing more and dimmer 
stars to be included in the observing program. 

The infrared detectors used in this study were 
gallium doped germanium bolometers operating at liquid 
helium temperatures. The earlier (1968) systems were 
pumped to just below the lambda point of liquid helium 
(2.2 °K) , whereas later improved dewar systems 
facilitated operation near 1.2 °K. 

The photometric systems of the 1968 vintage had 
only a single wavelength available. The effective 
wavelength was determined by an interference filter at 
liquid helium temperature. In these same dewars , the 
heat shield was cooled to approximately 150 °K by 
the escaping helium blow off gas. 

In addition to the 3.5y single wavelength systems, 
since 1970, it became possible to make observations at 
many wavelengths on the same night with the development 
of a multi-filter photometric system. The dewar utilized 
a ten position filter wheel in thermal contact with the 
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heat shield which was kept at about 80 °K by a liquid 
nitrogen container inside the dewar. In these systems, 
the basic set of standard filters available had 
effective wavelengths of 2.3, 3.5, 4.8, 8.6, 10.7, 

12. 2\x , and when atmospheric water vapor content 
allowed, 18y. 

As previously stated, as systems changed and 
improved, every effort was made to keep the inter- 
ference filters nearly identical, especially at 3.5y. 

Of course, this effort was not totally successful. 

The interference filters were purchased from Optical 
Coating Laboratories, Inc. (OCLI) and were originally 
intended to closely conform to Johnson's (1965) photo- 
metric systems. 

Equation 1 describes the method for determining 
the effective wavelength (Ae) of the photometric system. 

CO 

\ AT f (A) T a (A) D (a ) B ( A ) dA 

(U 

T p (A) T a (A) D ( A ) B (A ) dA 


T p (A) is the transmission of the interference filter 
obtained from a dual beam spectrophotometer. T.(A) 

li 

is the atmospheric transmission, D(A ) is the detector 
response and B (A ) is the energy distribution of the 


source . 
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In practice, the detector response D(A ) is 
assumed constant and independent of wavelength due 
to the nature of a bolometer. Also, B (A ) is 
approximately proportional to A ^ particularly long- 
ward of about 3y because the source energy distri- 
bution may be approximated by the Rayleigh-Jeans tail 
of a black body. The temperatures of most stellar 
energy sources exceed 2000 °K and their peak energy 
output would then be shortward of about 1 . 5y . 

Figure 2 is a dual beam spectrophotometer tracing 
of the transmission (T p ) of the interference filter 
used for most of the 3 . 5y measurements. This trans- 
mission curve was obtained with the filter at room 
temperature. When the interference filter is 
cooled to liquid nitrogen or liquid helium temperature, 
the bandpass has a tendency to shift toward shorter 
wavelengths. This shift may be as much as 0 . ly . 

The atmospheric infrared transmission (T A ) has 
been studied by Taylor and Jates (1957) . The trans- 
mission is a variable which depends on the amount of 
precipitable water vapor in the earth's atmosphere. 

In Minnesota, during the cold, dry winters, the water 
vapor content is typically less than 1 mm for which 
the infrared transmission is approximately constant 
between 3 and 4y . During the warm, humid summers, 
this HjO content may exceed 10 mm for which the 
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the transmission is greatly reduced shortward of 
about 3.3y. 

With allowances for these two aforementioned 
effects, the effective wavelength of the photometric 
system used for most of the observations was deter- 
mined to be 3.5y + 0.07 . 

A few observations at other wavelengths were 
obtained on some of the stars of this research 
program. Many different dewars , filters and photo- 
metric systems were employed but the systems had 
some common or equivalent filters such that the 
results should not be system dependent. 

The absolute calibration of the 3.5y filters and 
other wavelengths is based primarily on the work of 
Murdock (1972) , Johnson (1966) , Gehrz and Woolf (1971) , 
Gillett, Merrill and Stein (1971) and Gehrz (1972). 

The absolute calibrations of Gehrz and Woolf (1971) 
and Gillett et al. , (1971) have been used for the 3.5y 
effective wavelength in which zero magnitude corres— 
ponds to 7.1 x 10 watts/cm /y . Table I lists the 
effective wavelengths, Ae, and bandpass, AA, and the 
flux for zero magnitude for all of the filters used 
as a part of this research. 



WAVELENGTH (MICRONS) 



Figure 2. Dual beam spectrophotometer tracing of the 3.5 y filter. 
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Effective 


x e (u> 

2.3 

3.5 
4.8 

8.6 
10.7 
11.3 
12.2 
18 


Table I 

wavelength, bandpass, and absolute calibration 
of the Minnesota photometric system 


Flux for zero magnitude 
AAp F (A) in watts/cm 2 /y 


0.5 


1.1 


0.9 


0.9 


1.1 


2.0 


1.1 


4 


3.3 x 10 

7.1 x 10 

2.2 x 10 

2.3 x 10 
9.8 x 10 

8.0 x 10 

6.0 x 10 

1.3 x 10 


-14 

-15 

-15 

-16 

-17 

-17 

-17 

-17 
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IV 

SELECTION OF OBJECTS AND OBSERVATIONAL PROCEDURE 

From various sources, a list of the visually 
brightest semi-regular and Mira type variable stars 
was compiled. Preliminary measurements at 3.5y on 
some of these stars indicated that the proposed re- 
search project was indeed feasible, particularly for 
the cooler stars of later spectral type. From these 
survey observations and the sensitivity of the 
Minnesota infrared equipment available in 1968, the 
decision was made to observe only the brighter (at 
3.5y) variables. A semi-final list of program stars 
was then compiled which numbered about 35 stars 
distributed throughout the sky north of the celestial 
equator. The judgment to observe chiefly objects 
north of zero degrees declination was made to 
minimize the effects of atmospheric extinction. The 
program stars selected in this manner were mainly 
stars of spectral type M2 through M6 at maximum as 
well as several carbon stars and a few S stars. 

During the compilation of the program stars, it 
was not known by what magnitude, if any, nor on what 
time scale, other than the visual period, these red 
giants vary at infrared wavelengths. For this reason, 
it was decided to measure the selected objects as 
often as clear weather and available telescope time 
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permitted, hopefully every 5 to 10 days. Because Mira 
type variables have periods approximating one year, 
observations would be needed throughout the entire 
year, not just during the portion of the year when 
the star crosses the meridian between sunset and 
sunrise. In other words, daylight data were necessary 
although it had not been confirmed that accurate 
infrared observations could be obtained during daylight 
hours, especially when neither the object nor any field 
stars nearby could be visually identified. There was 
hope, however. 

Within the framework of the limitations previously 
described, the program stars were selected. The 
brighter infrared sources were chosen with aspirations 
of daytime observing but with the realization of the 
limitations imposed by the equipment. An effort was 
made to keep at least some representation of earlier 
(MO-M2 ) stars, and to maintain a somewhat uniform 
distribution of objects throughout the sky. However, 
stars near spectral type M6 at maximum tend to dominate 
the source list. 

The earliest observational data employing the 
inherently noisier detector systems was particularly 
difficult to obtain, because the observers were learning 
both infrared techniques and astronomy simultaneously. 

The original method of data acquisition consisted 
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of measuring a calibration star, a variable, the same 
calibration object, the same variable, etc. An 
evening's data were composed of observations of a few 
(-.5) calibration stars each measured 3 separate times 
and about five variable stars each observed at three 
distinct times. A typical integration time was 10 
minutes for each observation of each star while about 
10 minutes elapsed acquiring the next object. This 
procedure filled the entire observing run which lasted 
as long as clear skies allowed or. until the dewar 
exhausted its helium supply after about 8 hours at 
the telescope. 

The actual observing schedule and procedure evolved 
throughout the time span of data collection. By 1971, 
the procedure consisted of measuring each calibration 
star available once per session and each variable once 
per night; the implicit assumption being if the cali- 
brations were correct to a few percent then the vari- 
ables were accurate to that same value. When doubts 
existed, the specific stars were re-observed during 
the same session. The observational log consisted 
approximately of calibration, variable, new calibration 

star, new variable This method was employed 

so that any alterations in system sensitivity, atmos- 
pheric transmission, undetected clouds and other factors 
could hopefully be noted and, consequently, corrections 
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applied. 

About one year into the observing program, it 
was noted that Mira type stars do vary at 3 . 5y on a 
time scale compatible with the visual variations, 
although the amount of the change was considerably less, 
being about 1 magnitude at 3.5y compared to an average 
visual range of about 5 magnitudes. The original 
observing schedule of one data point per star every 
5-10 days was revised to one point every 10-20 days 
or, hopefully, about 30 points per cycle. 

The optimistic outlook was to follow the planned 
schedule for one or two cycles of each variable star. 
Such was not the case. The largest obstacle in the 
path of the proposed program was the lack of clear 
photometric quality weather in Minnesota. The second 
pitfall was the amount of available telescope time. 

The planned regular observations were purely at the 
mercy of the elements and nothing even approximating 
the planned data collection rate was achieved. Gaps 
in the observations would hopefully be filled by 
measurements during the next cycle. This latter 
assumption rested upon yet another hypothesis that 
the infrared magnitudes were reproducible from cycle 
to cycle. 

In an attempt to amass more complete data, daylight 
observations were attempted in the summer of 1969. By 
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this date, the O'Brien telescope had been debugged 
and tuned up to the extent that it could be accurately 
pointed to + 1 arc rain on the sky. These first day- 
time measurements were obtained with the 1-1/4 arc min 
beam diameter 3.5p system. After calibration, the 
telescope was set to the position of the variable to 
be measured. More often than not, the chart recorder 
indicated the star was within one of the two beams. 

If not, the telescope was driven in a raster pattern 
nearby until the signal was acquired from the star. 

The signal was maximized and the variable star measured. 
The chances of measuring a wrong star with the correct 
magnitude and position of the variable by this method 
were negligibly small. 

With this observational procedure, it was noted 
that daylight observations were only slightly more 
difficult than during the night. Drifts and variable 
offsets caused by solar heating of the telescope 
walls, defocusing of the mirror, and gradients in the 
sky background were minimal , at least at 3 . 5p with 
1-1/4 arc min diameter beams. In general, these 
daytime observations were limited to greater than 
about 20 degrees from the sun. 

Even with the improvements in detectors and 
systems since early 1968, some of the program stars 
were too dim at 3.5p to acquire during the sunlit 
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hours and gaps are obvious in the data. However, by 
1971, systems had progressed so that almost any of 
the program stars could be measured night or day by 
scanning the correct position and maximizing the 
3.5y signal. In fact, even with beams reduced to 
25 arc sec diameter , acquisition of an object was often 
faster using this scanning method than to visually 
identify the correct star field and then the proper 
variable. Although, when possible, visual verification 
was often attempted. 

Additional observations of some of the program 
stars were obtained at wavelengths other than 3 . 5y 
with some of the multi-filter systems. The basic 
procedure was the same as previously described for 
3.5y except that, due to the increased sky noise and 
the declining intensity of the source, the integration 

time necessary to achieve statistical significance was 

/ l 

increased atj the longer wavelengths . 

The multi-filter photometric system had an improved 
signal-to-noise ratio at 3.5y compared to the single 
wavelength dewar. However, because it employed a 
smaller beam diameter than the single 3. 5y system 
(25 arc sec vs. 75 arc sec) , the increased gain was 
partially offset by the poorer response profile and 
difficulties in accurate guiding. Daytime measurements 
were still feasible although the time needed to acquire 
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the stellar signal by scanning was increased. 

Data reduction was a relatively uncomplicated 
procedure. The digital output of the accumulator 
was recorded automatically. The object name, hour 
angle, universal time, and amplifier gain were 
recorded manually. The accumulator output was then 
returned to the laboratory where means and standard 
deviations of the mean could be computed. The total 
integration time on a specific star was quite 
subjective and depended on the signal-to-noise ratio, 
quality of the night, photometric system in use, 
infrared and visual seeing and transmission, as well 
as the philosophical bias of the observer. Generally, 
about 3 minutes of true integration time were deemed 
sufficient to reduce the statistical errors well 
below other uncertainties. 

The measured stellar fluxes were converted to 
magnitudes at unit air mass by equation 2. 

m, = m + e (AM, -AM )- 2.5 log (F,/F ) (2) 

1 S X S X 5 

Where and m g are the stellar magnitudes at 1.00 
air masses of the unknown and standard objects 
respectively ,e is the extinction coefficient in 
magnitudes per air mass (e<0) , AM, and AM are the air 
masses through which the unknown and standard are 
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observed while F. and F are the measured fluxes from 

1 S 

the unknown and the standard. 

The air mass through which a star was observed 
was calculated from the standard formula 3 

AM = sec Z (3) 

where Z is the zenith angle. Because the program stars 
were always at zenith angles less than 60° and generally 
less than 45° equation (3) is an adequate approximation. 

The atmospheric extinction in the 3 to 4y window 
was measured during each observing session whenever 
possible. The mean value during the four years of 
data collection was 0.2 magnitudes per air mass at 
the 3.5y effective wavelength. After the first two 
years of "observation , this value was checked less 
frequently but was found to be the same within 
observational errors. 

This value of 0.2 magnitudes per air mass for 
the atmospheric extinction at 3.5y was adopted and used 
for all of the extinction corrections on all of the 
data. Due to the fact that the program stars generally 
had positive declinations and were measured as near 
the meridian as possible, for the 45° latitude of the 
Minnesota telescope, the air masses were always less 
than 1.5 and the extinction correction was less than 
0.10 mag. The 30% variations of e, as measured from 
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night to night and season to season, therefore, are 
able to produce only a + 0.03 mag .uncertainty in the 
final results. 

From the preliminary and first year's measure- 
ments , it was pleasantly discovered that observations 
of brighter stars on a given night were reproducible 
to about 2%. These errors are assumed to be caused 
by inaccurate guiding, beam positioning, telescope 
tracking, fluctuations in atmospheric transmission, 
poor seeing, non-linearity of the photometric system, 
statistical fluctuations, and other assorted gremlins. 
For sources whose 3 . 5y magnitude was dimmer than about 
+1™5 , the chief uncertainty was statistical, due to 
the poorer signal-to-noise ratio and limited inte- 
gration time. For stars below about 1.5 air masses, 
the extinction must be more accurately determined 
for each night than was attempted. 

A separate discrepancy in the absolute flux 
calculated for the variables may result if the assumed 
flux for zero magnitude is not correct. 

The zero point of the Minnesota photometric system 
was defined by the assumption that a Lyra (AOV) was 
zero magnitude at all wavelengths. The magnitudes of 
the calibration stars for the preliminary observations 
were adopted as published by Johnson (1964). At 
that time, it was discovered that the magnitudes 
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calculated for some of the calibration stars did not 
agree to within + 10% of the quoted values. In 
general, stars of later spectral type were brighter 
than as measured by the Arizona equipment, possibly 
indicating the difference in effective wavelength as 
defined by the two groups. Due to this possible 
systematic difference and, in order to define more 
calibration stars, a new list of standards was compiled. 
This new list was comprised of stars of later spectral 
type, bright at 3.5y, and whose positions on the 
celestial sphere were more compatible to the distribution 
of variables in the observational program. This new 
calibration list was based on a Lyra still defined as 
zero magnitude, and the standards were measured 
directly with respect to Vega if possible. 

Table II lists the primary standard stars by 
their HR number, name, spectral type and 3 . 5y 
magnitude defined by the Minnesota system. Table III 
compiles the same quantities for the secondary 
calibration stars. All of the magnitudes listed 
include corrections for atmospheric extinction and are 
referenced to unit air masses. 

When the results of Minnesota photometry can be 
compared with magnitudes published by other groups 
(Lee, 1970; Gillett, Merrill and Stein, 1971; for 
example) there exists surprisingly good agreement, 
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3.5y 

Magnitudes 

Table II 

of Primary Standard 

Stars 

HR 

Name 

Spectral 
Type * 3 

.5u Magnitude 

337 

8 And 

MO III 

-2.05 

1457 

a Tau 

K5III 

-3.02 

1708 

a Aur 

G8III7+F 

-1.92 

2943 

a CMi 

F5IV 

-0.76 

2990 

8 Gem 

KO III 

-1.23 

4058 

2 _ 

Y Leo 

G7III 

-0.79 

4301 

a UMa 

KO II-III 

-0.77 

5340 

a Boo 

K2IIIp 

-3.12 

6705 

y Dra 

K5III 

-1.44 

7001 

a Lyra 

AOV 

0.00 

8775 

8 Peg 

M2II-III 

-2.45 


* Hoffleit, 1964 



Table III 
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3.5\i Magnitudes of Secondary Standard Stars 


Spectral 

HR Name Type* 3 . 5u Magnitude 


168 

a 

Cas 

KO 

II-III 

-0.40 

603 

Y' 

' And 

K3 

II 

-1.00 

617 

a 

Ari 

K2 

III 

-0.77 

911 

a 

Cet 

M2 

III 

-1.90 

921 

P 

Per 

M4 

II-III 

-2.20 

1577 

l 

Aur 

K3 

II 

-0.85 

2061 

a 

Ori 

M2 

lab 

-4.54 

2091 

TT 

Aur 

M3 

II 

1 

• 

O 

2216 

n 

Gem 

M3 

III 

-1.64 

2286 

y 

Gem 

M3 

III 

-2.10 

3705 

a 

Lyu 

MO 

III 

-0.76 

4069 

y 

UMa 

MO 

III 

O 

0 

* 

.H 

1 

4434 

A 

Dra 

MO 

III 

-0.28 

4910 

6 

Vir 

M3 

III 

-1.40 

5506 

£ 

Boo 

KO 

II-III 

-0.07 

6406 

a 1 

Her 

M5 

II 

-3.70 

7139 

6 2 

> 

Lyra 

M4 

II 

-1.41 

7157 

R 

Lyra 

M5 

III 

-2.28 

7417 

P 

Cyg 

K5II7+B? 

0.05 

7525 

Y 

Aql 

K3 

II 

-0.75 

7536 

5 

Sge 

M2II7+B? 

-0.96 

7949 

£ 

Cyg 

KO 

III 

0.03 

8308 

£ 

Peg 

K2 

lb 

00 

o 

8316 

y 

Cep 

M2 

la 

-2.23 

9064 

if 

Peg 

M3 

III 

r* 

« — i 

o 


* Hoffleit, 1964 



generally to within 2-3%, even though effective 
wavelengths may differ somewhat. 
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V 

OBSERVATIONAL RESULTS AT 3.5y 

Table 4 summarizes the characteristics of the 
variable stars of this research project for which more 
than a few (~5) reliable observations have been secured. 
Column 1 lists the star name by constellation, while 
columns 2 and 3 tabulate the galactic coordinates of 
longitude and latitude, respectively (Kukarkin, et al., 
1958) . Column 4 quotes the type of variation as Mira 
type (M) or semi-regular (SR) (Kukarkin, et al. , 1958). 
Columns 5 and 6 list the average periods cited by 
Kukarkin, et al. , (1958) and Campbell (1955), whereas 
column 7 lists the adopted period used for the phase 
calculations. Column 8 lists the mean spectral types 
at maximum (Keenan, 1966) or the spectral type 
(Kukarkin, et al., 1958) and column 9 tabulates the 
minimum spectral types of Lockwood and Wing (1971) or 
Kukarkin, et al., (1958). Columns 10, 11 and 12 quote 
the maximum and minimum mean visual magnitude and 
range determined by Campbell (1955) . Last, but not 
least, columns 13, 14 and 15 tabulate the observed 
3 • 5y magnitude extrema and range as determined by 
this present project. The extrema of the last columns 
represent the brightest and dimmest reliable measure- 
ments at 3.5y and probably are not a true measure of the 
mean behavior of the stellar fluctuations. When the 



Table IV 

Characteristics of the Program Stars 





Type of 

Period 

Adopted 

Spectral 

Spectral 

Mean 

Visual 

+ 


3 . 5 j 



l 

b 

Variation 

Days 


Period 

Type at 

Type at 

Magnitude 



Magnitude 


Name 

* 

# 

* 

A 

+ 

Days 

Maximum 

Minimum 

Max 

Min 

Range 

Max 

Min 

Range 

R And 

085° 

-24 0 

M^Mira 

409.4 

409.0 

410 

54 , 7e + 

- 

6?9 

14T‘3 

7?4 

-1?3 

-0^1 

1?2 

W And 

107 

-15 

M 

397.0 

397.0 

392 

S8,2e + 

M8e * 

7.4 

13.7 

6.3 

-0.6 

0.0 

0,6 

R Aql 

010 

-01 

M 

300.3 

306.2 

283 

M6,5e + 

M8.5 : 

6.1 

11.5 

5.4 

-1.5 

-1.1 

0.4 

R Aqr 

037 

-71 

H 

386.9 

386.4 

387 

M7e+pec * 

- 

6.5 

10.3 

3.8 

-2.3 

(-1.9) 

- 

P Aur 

124 

+ 10 

M 

458.9 

461.8 

454 

M6.5e + 

M9.0 : 

7. 7 

13.5 

5.8 

(-1.6) 

-1.1 

* 

R Boo 

005 

65 

M 

223.3 

223.4 

224 

K4.5e + 

M8.0 : 

7.2 

12.3 

5.1 

1.0 

1.7 

0.7 

R Cnc 

180 

26 

M 

362.1 

2 62.2 

360 

M6.5e + 

M8e * 

6.8 

11.2 

4.4 

-1.3 

(-1.0) 

- 

T Cnc 

175 

38 

SRa 

482.4 

- 

- 

K3(C4 5 ) * 

- 

7.6 

10.5 

* 

0.2 

0.4 

(0.2) 

R Cas 

083 

-10 

M 

431.2 

429.6 

422 

M7e + 

M10 s 

7.0 

12.6 

5.6 

(-2.6) 

-2.1 

- 

T Cas 

08? 

-07 

M 

445.0 

445.0 

445 

M7.5e + 

M8.5 j 

7.9 

11.9 

4.0 

(-1.8) 

-1.4 

- 

WZ Cas 

034 

-02 

SRa 

1S6.0 

- 

- 

Kp(C9 L )* 

- 

9 * 4p 

11. Ip 

* 

-0.1 

0.4 

0.5 

S Cep 

081 

19 

M 

437.5 

484. S 

482 

NBe(C7 4 e) * 

- 

8.3 

11.0 

2.7 

-1.8 

-1.2 

0.6 

T Cep 

072 

14 

M 

389.7 

392.9 

382 

M6.5e + 

’48 . 5 s 

6.0 

10.3 

4.3 

-2.4 

-2.2 

(0 . 2 i 

0 Cet 

137 

-57 

M 

331.6 

331.4 

332 

M5.5e + 

M9.0 : 

3.5 

9.2 

5.7 

-3.7 

-2.9 

0.6 

S CrB 

016 

56 

M 

360.7 

359.5 

359 

M6.5e + 

M8e * 

7.3 

12. S 

5.6 

-0.8 

-0.3 

0.5 

V CrB 

030 

50 

M 

357.8 

357.6 

358 

N2(C6 2 e) * 

- 

7.5 

11.0 

3.5 

0.3 

1.1 

0.6 

R Cyg 

050 

13 

M 

426.3 

426.3 

427 

S3,9e + 

S6 t 8e * 

7.5 

13.9 

6.4 

0.2 

1.4 

1.2 






Type of 

Period 



l 

b 

Variation 

Days 


Name 

* 

a 

* 



U Cyg 

052 

06 

M 

464.7 

462.9 

RS Cyg 

044 

02 

SRa 

417.8 

419.0 

kt cyg 

049 

11 

M 

190.4 

190.7 

X Cyg 

036 

02 

M 

406.9 

409.0 

RY Dra 

088 

51 

SRb 

172.5 

- 

U Her 

003 

39 

M 

406.0 

405.3 

R Leo 

192 

46 

M 

312.6 

312.8 

R LMi 

133 

51 

M 

372.3 

372.2 

R Lap 

182 

-30 

M 

432.5 

- 

x oph 

007 

06 

M 

334.2 

334.3 

U Ori 

156 

-01 

M 

372.3 

373.2 

w ori 

167 

-21 

SRb 

212 

- 

R Ser 

354 

45 

M 

356.8 

356.6 

R Tri 

115 

-23 

K 

266.4 

266.0 

R UMa 

105 

45 

M 

301.7 

301.0 

R Vir 

265 

69 

M 

145.6 

145.4 

SS Vir 

259 

63 

M 

354.7 

355.9 

~*g»lraT!ri^ 

et al., 

1958. 


♦Keenan, 1966. 


Table IV (continued) 


Adopted 

Spectral 

Spectral 

Mean Visual 

■f 


3- Sii 


Period 

Type at 

Type of 


Magnitude 


Magnitude 


Days 

Maximum 

Minimum 

Max 

Min 

Range 

Max 

Min 

Rar.ee 

455 

Npe(C7 2 e) * 

C9 2 * 

7-2 

10.7 

3.5 

-0.1 

1.0 

1.1 

418 

NOpe (CS^e) * 

- 

7.2 

9.0 

1.8 

0. 0 

1.0 

(0.2) 

190 

M2 . 5e + 

M8.2 i 

7.3 

11.8 

4.5 

2.5 

3.0 

(c.s: 

411 

S7,2e + 

M9.0 i 

5.2 

13.4 

8.2 

-3.1 

-2.4 

0.7 

- 

tf4p(C3 4 ) * 

- 

9 . 4p 

11. 4p 

» 

-0.6 

-0.4 

(0.2) 

413 

M6.5e ♦ 

M9.0 : 

7.5 

12.5 

5.0 

(-1.1) 

-0.6 

- 

310 

M7e + 

M9.0 : 

5.9 

10.1 

4.2 

(-3.4) 

-2.9 

- 

372 

M7e + 

M8.8 ; 

7.1 

12.7 

5.6 

-1.4 

(-1.0) 

- 

- 

N6e(C7 2 e) * 

- 

5.9 

10.5 

4.6 

- 

- 

- 

332 

M6.5e + 

M8.5 : 

6.8 

8.8 

2.0 

-1.6 

-1.2 

0.4 

372 

M6.5e + 

M9.0 : 

6.3 

12.0 

5.7 

-1.7 

-1.0 

0.7 

- 

N5(C5 3 ) * 

- 

5.9 

7.7 

* 

-1.0 

-0.8 

(0.2) 

352 

MG . 5e + 

M8e * 

6.9 

13.4 

6.5 

(-0.4) 

0.0 

- 

267 

M4e + 

M8.2 s 

6.2 

11.7 

5.5 

0.1 

0.6 

0.5 

302 

M4.5e + 

M6e * 

7.5 

13.0 

5.5 

(0.5) 

1.0 

w 

146 

M4.5e + 

MSe * 

6.9 

11.5 

4.6 

1.4 

(1.6) 

- 

355 

Ne(C6 3 e) * 

_ 

6.8 

8.9 

2.1 

0.0 

(0.3) 

- 


: Lockwood and Wing,. 1971, 


+ Campbell (1955) 



maxima or minima are enclosed in parentheses, it 
indicates that the observed extrema are caused by 
limitations in the available data and therefore do 
not depict the true stellar range of brightness. 

When the 3 . 5y range is on the order of the un- 
certainties of measurement (~0 . 2 ) , it is also enclosed 
in parentheses and indicates that the implied altera- 
tions may not be real. 

The adopted periods of column 7 , as used in the 
phase calculations, may differ somewhat from the aver- 
age periods of Campbell (1955) in column 6. The 
adopted periods were calculated from visual data 
kindly supplied by Mayall (1972) from Julian dates 
2439000 through 2441500. These periods represent an 
average period of the star over the inclusive dates 
and are intended to depict the present, rather than 
the past, average behavior. Only in one notable 
exception do the adopted periods differ from Campbell' 
averages by more than 11 days or more than about 3%. 
Generally, the difference in the two periods is less 
than about 5 days and therefore represents an error of 
less than 2%. Differences of 0.01 in the phase 
calculations are much less than the cycle to cycle 
variations in the period. The notable exception is 
R Aql , which is known to have a decreasing period as 
previously stated. 
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Appendix I summarizes all of the program stars' 
measurements as observed with the Minnesota photo- 
metric system's 3.5y effective wavelength. The 
headings for each star include the star name and the 
average period used to determine the phase. Column 1 
lists the Julian date (rounded to the nearest whole 
day) of the observation, column 2 represents the phase 
with respect to the visual maximum, column 3 tabulates 
the magnitude at 3 . 5y and column 4 contains an 
estimate of the quality of the data. 

Appendix II tabulates the 3.5u observations for 
sources other than the final program stars which have 
been measured more than one time. It is included 
only as a summary. Column 1 lists the star name 
(alphabetically by constellation); column 2, the 
Julian date of observation; columns 3 and 4, the 3.5p 
magnitude and quality estimate, similar to appendix I; 
columns 5, 6, 7, 8 list the type of variability, 
period, spectral type and visual range, respectively. 
The period, spectral type and visual magnitudes are 
from Campbell (1955) , if available, and from 
Kukarkin, et al. (1958) , when noted. 

The phase of the star at the time of the infrared 
observation is measured with respect to the time of 
visual maximum. As a general rule, the phase is 
measured from zero at one visual maximum through 
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phase + 1.00 at the ensuing maximum, thus making 
the phase always positive and ranging between 0.00 
and +1.00. However, cycle to cycle differences exist 
in both the visual and infrared observed magnitudes. 

A visual maximum may be brighter than average, and, 
if so, there exists a tendency for this particular 
extremum to be earlier than average (Harrington, 

1965) implying a shorter observed period for that 
cycle. Phase calculations from 0 through +1 tend 
to split this brighter maximum into two different 
cycles. In order to reduce this apparent discontinuity 
at zero phase, it was decided to measure phases on 
either side of the visual maximum, from phase -0.50 
to 0.00 to +0.50 and then recycle. This method of 
computation stresses the time of visual maximum and 
confines the unusual brightness to only one cycle and 
not over two cycles. 

The 3 . 5y magnitudes tabulated in column 3 
include the effects of atmospheric extinction and 
refer to unit air mass. 

The numbers listed in column 4 (Quality) are 
a measure of the performance of the photometric 
system's noise, beam size, and stability throughout 
the observing session. These numbers also represent 
constancy of the atmospheric transmission, cloud cover, 
if any, visual and infrared seeing, as well as assorted 
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other objective and subjective decisions. The major 
factor in determining the quality of the observations 
was the accuracy to which the calibrations stars 
were measurable and reproducible during the observing 
run . 

On an excellent photometric night, the atmos- 
pheric extinction, as measured, was found to agree 
within about 20% of the adopted 0.2 magnitudes per 
air mass. Also, there were few or no clouds present, 
and the detector system was stable to within a few 
(~2) percent throughout the observing session. On 
these excellent nights, the calibration stars were 
reproducible to about 2 or 3 percent, and the night 
was rated as quality 5. 

Quality 4 implies a good photometric night. 

Some of the aforementioned criteria may have influenced 
the results. For example, the calibration stars were 
accurate to 2 or 3 percent after corrections were 
applied for changing sensitivity or transmission. 
Quality 3 implies fair reproducibility, while 2 
represents a poor night. Quality 1 data are not 
included because, for a multitude of reasons, the 
results are not reliable to within an estimated 50%. 

It is appealing to assign overly optimistic 
errors to the excellent nights. During these observing 
sessions, as previously stated, the calibration stars 
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were accurate to about 2%. It is tempting to assume 
that the variable stars are accurate to the same 2%. 
Even though this may be often a reality, a conservative 
estimate would be doubling this 2%; therefore, nights 
of quality 5 were given 4% uncertainties. 

Along a similar vein, quality estimate 4 nights 
were most likely accurate to better than 5% on the 
calibration stars after the various correction 
factors were applied to the data. However, the 
errors are conservatively estimated to be between 
5 and 10%. An arbitrary error of 10-20% is assigned 
to quality 3 nights due to the presence of clouds , 
equipment problems, and more or larger applied 
correction factors. Last, but not least, quality 
estimate 2 represents basically poor data reliability 
of 20 to 50% errors. 

At this point, a brief discussion of the sources 
of errors and uncertainties involved in limiting the 
accuracy of the 3.5u magnitudes is in order. The 
chief difficulty in obtaining observational data, at 
least in Minnesota, is the lack of clear or partly 
clear weather — especially during the daylight hours. 
One major producer of poor data quality is unseen, 
thin cirrus cloud cover. At the 3 . 5y wavelength, 
these clouds do not introduce additional noise; they 
simply reduce the signal. In fact, thin uniform 



42 


overcast cirrus may reduce the observed signal level 
by 20 to 30% and this reduction may go unnoticed for 
some time. If this reduction did occur during an 
observing session, it was usually noticeable upon 
final analysis of the data and this correspondingly 
poor data was omitted. Clouds, fog, and haze are 
classified as major sources of observational uncertainty 
Many possible causes of minor errors do also 
exist. Minor errors are the 1 or 2 percent variety 
which may go undetected for some time, quite possibly 
forever. One gremlin is the changing sensitivity of 
the bolometer. This changing sensitivity is due to 
oscillations in the vacuum system and/or a long time 
constant for the dewar to reach equilibrium. Additional 
uncertainties may be caused by non-linearity in the 
electronics, fogging over of a dewar window in humid 
weather, or cracking of the window in very cold 
conditions. More sources of minor error include 
positioning of the telescope to obtain maximum source 
signal strength, non-uniformity of response over the 
area of the detector, and non-linearity of the detector 
for very bright sources. As previously mentioned, the 
changing atmospheric transmission and the associated 
extinction do also affect the accuracy of the results. 
Last, but not least, are the Minnesota winters where 
the temperature does occasionally dip to -40 F. This 
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may cause severe observational problems; however, the 
cold weather observing technique is an entirely 
separate tale to be told to one's grandchildren. 

Many of the possible sources of error can be 
eliminated or at least minimized. The ones which 
cannot be eliminated or minimized may hopefully be 
accounted for and corrections applied. As also 
stated previously, the final errors were estimated 
from the quality and degree of reproducibility of the 
data and not from the possible problems. 

Statistical errors are the major producer of 
uncertainty for objects with a low signal-to-noise 
ratio and a limited integration time. For the 
earlier (1968-1970) observations, the minimum 
detectable signal, in about twenty minutes of real 
integration time, was a 3.5y magnitude of about +4 ra 
at the three standard deviation limit. Since 1971, 
this limit has improved to about +7 m in the same 
integration time. Most of the program stars are 
brighter than +l m and the statistical errors are not 
the prime source of uncertainty. But for objects near 
2 m , the statistics start to become important, parti- 
cularly for the daylight observations; therefore, for 
objects with a 3.5y magnitude near +2 m , the assigned 
errors are almost always greater than 20%. 

Figure 3 plots the visual variation and the 3.5y 
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J.D. 2440000 + 

Figure 3a. Continuous visual and 3.5y light curves for the mild S star x Cyg(S7,2e). 

Typical errors are +5% for the 3.5y points. The arrows denote the times 
of visual extrema. 
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variability as functions of Julian date for x Cyg, 

R Cyg, R Aql , and U Cyg. These stars were selected 
because their data are among the most complete and 
continuous. The visual curves represent an eye 
estimate of the raw observational data supplied by 
the AAVSO. Representative error bars for the 3.5p 
observations are typically about + 0™0 5, unless 
otherwise noted. The times of visual maximum 
(AAVSO) and minimum are noted by arrows on the graphs. 

Figure 4 graphs the 3 . 5p magnitude as a function 
of phase for all of the M stars in the program for 
which phases could be computed. The distinct light 
cycles are coded as follows: + = cycle 1, o = 2, 

• = 3, x - 4, A = 5, □= 6. Figure 4a plots those 
stars with periods less than 300 days, figure 4b does 
the same for stars whose periods lie between 301 and 
350 days; figure 4c is the same except for periods 
between 351 and 400 days and figure 4d includes stars 
whose periods exceed 400 days. Figure 5 is similar to 
figure 4 but includes only S type stars and figure 6 
graphs the data for carbon stars. 

It is obvious that some stars in the original 
program listed in table IV are not plotted in figures 
4, 5 or 6. These omissions of T Cnc , WZ Cas , RT Cyg, 

RY Dra, and W Ori merit further comment: 

rl 

T Cnc [N3 (C4,. ) ,482 , IRC+20 207 ] is a semi-regular variable 
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Figure 4a. 3 . 5y light curves as a function of visual 

phase for the M stars with periods less than 
300 days. The distinct cycles are coded as 
follows: +=cycle 1 , 0=2, t=3, x=4 , A=5,Q=6. 
The errors in the measurements are estimated 
to be between 5 and 10%. The magnitudes en- 
closed in circles are of poorer quality. The 
curve drawn through the measurements repre- 
sents a best "eye" estimate as to the shape 
and mean value of the 3.5y variations. 



MAGNITUDE 


50 
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Figure 4b. Similar to figure 4a but for M stars with 
periods between 301 and 350 days. 
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Figure 4c. Similar to figure 4a but for M stars with 
periods between 351 and 400 days. 
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Figure 4d. Similar to figure 4a but for M stars 
with periods greater than 400 days. 
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Figure 5. Similar to figure 4a but for S stars. 
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Figure 6. Similar to figure 4a but for carbon stars. 
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whose visual data from the AAVSO were deemed in- 
sufficient to produce Julian dates of visual maxima. 
The scarcity of reliable 3.5y data points, which seem 
to exhibit little if any variation other than 
observational errors, allows little to be stated about 
variability and phase. The adopted 3 . 5y magnitude is 
+ 0^26 + 0^02 (1 standard deviation of the mean) from 
12 data points. 

WZ Cas [Np (C9-^) , 186 d , IRC+60433] is another 
SRa variable. Again, the AAVSO visual data were 
considered insufficient for determining dates of 
visual maxima. However, in contrast to other 
omitted stars, WZ Cas does seem to exhibit a real 
variation at the 3 . 5y wavelength. 

RT Cyg [M2.5e, 191 d , IRC+50306] , in contrast 
to the previous stars, is a Mira type variable. The 
AAVSO data well demonstrated the visual variations, 
but at 3 . 5y this star was too dim for reliable data. 

The adopted 3 . 5y magnitude is +2™1 A + 0^07 from 9 
observations . 

RY Dra [N4p (C3 4 ) , 172 d 5, IRC+70116 , SRb] : 

No visual observations from the AAVSO were available. 
The 3.5y magnitudes are scarce and range between -0 1 ?40 
to -0™59 for good quality observations. The adopted 
mean 3.5y magnitude is -0^52 + 0^02 . 

W Ori [N5 (C5 3 ) , 212 d , IRC+00066, SRb]: The AAVSO 
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visual data were again inadequate for days of visual 
maxima. The reliable 3 . 5u data were sparse and 
ranged between -0?82 to -l^Ol. The adopted mean 
magnitude was determined to be -0™94 + 0™02. 

Figure 7 is a plot of the 3.5u observational 
magnitudes as a function of visual phase for x Cyg f 
the star with the most complete and continuous data. 

This is the same data as in figure 5 but on an 
expanded magnitude and phase scale. From this curve, 
the reproducibility and quality of the data and the 
cycle to cycle changes in magnitude and shape are 
more easily seen. 

The solid lines drawn on the preceding figures 
represent a best eye estimate to the data. The 
attempt was made to determine the shape of the 3.5y 
light curve from observations on the same cycle. This 
general shape for each stellar curve was then adjusted 
(again, by eye) to represent a mean magnitude averaged 
over the separate cycles. Of course, data of poor 
quality were given less weight. The solid curves 
thus depict the best eye estimate to the available 
data, and represent shape and mean magnitude as functions 
of phase. It was felt that a more sophisticated^ 
analysis of the data, such as a least squares method, 
would not produce more meaningful curves because of 
the apparent cycle to cycle changes in magnitude 



- 0.5 - 0.4 - 0.3 - 0.2 - 0.1 0.0 0.1 0.2 0.3 0.4 0.5 


Figure 7. 
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3.5y light curve versus visual phase for the mild S star 
v Cvg(S7,2e) on an expanded magnitude and phase scale. 





Table V 

Adopted 3.5y Magnitudes as a Function of Phase 


Phase 

R And 

W And 

R Aql 

R Aqr 

R Aur 

R BOO 

R Cnc 

R Cas 

T Cas 

-0.50 

-0.72 

-0.18 

-1.16 

- 

-1.11 

1.52 

- 

-2.27 

-1.50 

. 45 

-0.58 

-0.12 

-1.25 

- 

-1.07 

1.57 

- 

-2.19 

-1.54 

.40 

-0.46 

-0.08 

-1.29 

- 

-1.07 

1.59 

- 

-2.15 

-1.55 

.35 

-0.36 

-0.04 

-1.30 

- 

-1.11 

1.60 

- 

-2.16 

-1.54 

.30 

-0.29 

-0.03 

-1.30 

- 

-1.13 

1.60 

- 

-2.20 

-1.53 

.25 

-0.27 

- 

-1.28 

- 

-1.13 

1.59 

- 

-2.25 

-1.50 

.20 

-0.28 

- 

-1.24 

-1.9 4 

-1.13 

1.58 

- 

-2.30 

-1.46 

.15 

-0.32 

- 

-1.22 

-1.97 

-1.14 

1.53 

-1.0 3 

-2.35 

-1.41 

.10 

-0.41 

- 

-1.25 

-2.03 

-1.17 

1.46 

-1.10 

-2.39 

-1.41 

-0.05 

-0.60 

- 

-1.29 

-2.13 

-1.25 

1.39 

-1.18 

-2.44 

-1.45 

0.00 

- 

-0.49 

-1.34 

-2.22 

-1.38 

1.31 

-1.25 

-2.51 

-1.54 

+0.05 

- 

-0.55 

-1.40 

-2.28 

-1.49 

1.24 

-1.29 

-2.59 

-1.65 

.10 

- 

-0.60 

-1.46 

-2.30 

-1.56 

1.19 

-1.32 

-2.66 

-1.72 

.15 

- 

-0.61 

-1.50 

-2.29 

-1.58 

1.17 

-1.33 

- 

-1.78 

.20 

-1.30 

-0.59 

-1.50 

-2.23 

-1.56 

1.19 

-1.32 

- 

-1.77 

.25 

-1.25 

-0.55 

-1.46 

-2.15 

-1.51 

1.24 

-1.31 

-2.60 

- 
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Table V (continued) 


Phase 

R And 

W And 

R Aql 

R Aqr 

R Aur 

R Boo 

R Cnc 

R Cas 

T Cas 

.30 

-1.17 

-0.50 

-1.40 

- 

-1.45 

1.29 

-1.26 

-2.53 

- 

.35 

-1.08 

-0.45 

-1.32 

- 

-1.38 

1.35 

-1.20 

-2.46 

-1.61 

.40 

-0.99 

-0.37 

-1.23 

- 

-1.30 

1.40 

-1.11 

-2.40 

-1.55 

. 45 

-0.86 

-0.27 

-1.15 

- 

-1.20 

1.46 

-1.00 

-2.33 

-1.49 

+ 0.50 

-0.72 

-0.18 

-1.16 

- 

-1.11 

1.52 

- 

-2.27 

-1.50 


Phase 

S Cep 

T Cep 

0 Cet 

S CrB 

V CrB 

R Cyg 

U Cyg 

-0.50 

-1.20 

-2.35 

-3.06 

-0.38 

0.96 

0.90 

1.01 

.45 

-1.23 

-2.30 

-3.05 

-0.32 

1.01 

1.10 

1.04 

.40 

-1.25 

-2.26 

-3.07 

-0.31 

1.05 

1.25 

0.92 

.35 

-1.25 

-2.24 

-3.11 

-0.33 

1.07 

1.37 

0.74 

. 30 

-1.27 

-2.22 

-3.13 

-0.36 

1.05 

1.41 

0.65 

.25 

-1.33 

-2.21 

-3.16 

-0.40 

1.00 

1.40 

0.61 

.20 

-1.41 

-2.20 

-3.20 

-0.43 

0.85 

1.34 

0.58 

.15 

-1.51 

-2.20 

-3.24 

-0.45 

0 . 60 

1.20 

0.54 

.10 

-1.62 

-2.20 

- * -3.29 

-0.47 

0.43 

0.86 

0.41 

-0.05 

-1.71 

-2.20 

-3.33 

-0.50 

- 

0.62 

0.22 

0.00 

-1.79 

-2.21 

-3.39 

-0.58 

- 

0.44 

0.02 
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Table V (continued) 


Phase 

S Cep 

T Cep 

0 Cet 

S CrB 

V CrB 

R cyg 

U Cyg 

+0.05 

-1.81 

-2.23 

-3.46 

-0.65 

- 

0.30 

-0.06 

.10 

-1.79 

-2.25 

-3.51 

-0.73 

- 

0.23 

-0.08 

.15 

-1.63 

-2.27 

-3.54 

-0.79 

0.33 

0.20 

-0.05 

.20 

-1.53 

-2.29 

-3.54 

-0.80 

0.41 

0.21 

0.04 

.25 

-1.44 

-2.31 

-3.52 

-0.78 

0.51 

0.26 

0.16 

.30 

-1.33 

-2.33 

-3.47 

-0.73 

0.60 

0.34 

0.26 

.35 

-1.21 

-2.35 

-3.38 

-0.66 

0.70 

0.43 

0.39 

.40 

-1.17 

-2.37 

-3.25 

-0.58 

0.79 

0.57 

0.58 

.45 

-1.17 

-2.37 

-3.13 

-0.48 

0.88 

0.71 

0.80 

+0.50 

-1.20 

-2.35 

-3.06 

-0.38 

0.96 

0.90 

1.01 


Phase 

X Cyg 

U Her 

R Leo 

R LMi 

X Oph 

U Ori 

R Ser 

R Tri 

-0.50 

-2.50 

-0.69 

-3.09 

- 

-1.20 

-1.16 

-0.08 

0.38 

.45 

-2 . 41 

-0.69 

-3.03 

- 

-1.20 

-1.08 

-0.07 

0.42 

. 40 

-2 . 39 

-0.73 

-3.03 

- 

-1.22 

-1.00 

-0.08 

0.45 

.35 

-2.45 

-0.80 

-3.06 

- 

-1.25 

-0.98 

-0.10 

0.45 

.30 

-2.53 

-0.86 

-3.08 

- 

-1.27 

-0.99 

-0.13 

0.44 
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Table V (continued) 


Phase 

X Cyg 

U Her 

R Leo 

.25 

-2.55 

-0.89 

-3.10 

.20 

-2.53 

-0.90 

-3.11 

.15 

-2.50 

-0.92 

-3.13 

.10 

-2.55 

-0.97 

-3.16 

in 

o 

o 

1 

-2 . 70 

-1.04 

-3.19 

0.00 

-2.90 

-1.13 

-3.24 

+ 0.05 

-3.01 

- 

-3.29 

.10 

-3.09 

- 

-3.35 

.15 

-3.11 

- 

-3.41 

.20 

-3.11 

-1.09 

-3.44 

.25 

-3.07 

-1.03 

-3.45 

.30 

-2.98 

-0.97 

-3.44 

. 35 

- 2.87 

-0.90 

-3.41 

.40 

- 2.73 

-0.84 

-3.35 

.45 

-2 . 61 

-0.76 

-3.26 

+ 0.50 

-2.50 

-0.69 

-3.09 


R LMi 

X Oph 

U Ori 

R Ser 

R Tri 

-1.08 

-1.30 

-1.02 

-0.15 

0.44 

1 

• 

O 

-1,33 

-1.07 

-0.19 

0.44 

-1.09 

-1.35 

-1.15 

-0.22 

0.45 

-1.12 

-1.36 

-1.24 

-0.26 

0.44 

-1.17 

-1.38 

-1.33 

-0.30 

0.41 

-1.22 

-1,42 

-1.44 

-0.34 

0.34 

-1.29 

-1.46 

-1.55 

- 

0.28 

-1.37 

-1.50 

-1.70 

- 

0.21 

-1.40 

-1.52 

-1.77 

- 

0.15 

l 

* 

o 

-1.54 

-1.76 

- 

0.14 

-1.38 

-1.51 

-1.65 

-0.44 

0.18 

-1.30 

-1.46 

-1.51 

-0.38 

0.23 

-1.22 

-1.41 

-1.40 

-0.31 

0.27 

-1.13 

-1.34 

-1.32 

-0.24 

0.31 

- 

-1.26 

-1.24 

-0.16 

0.35 

— 

-1.20 

-1.16 

-0.08 

0 .38 



Table V (continued) 


Phase 

R UMa 

R Vir 

SS V; 

-0.50 

0.95 

1.65 

- 

.45 

1.01 

- 

- 

.40 

1.01 

- 

- 

.35 

0.94 

- 

- 

.30 

0.88 

- 

- 

.25 

0.82 

- 

- 

.20 

0.76 

- 

- 

.15 

0.71 

- 

- 

.10 

0.65 

- 

0.30 

-0.05 

0.60 

- 

0.26 

0.00 

0.55 

1.58 

0.20 

+ 0.05 

0.50 

1.50 

0.13 

.10 

0.47 

1.43 

0.08 

.15 

0.45 

1.39 

0.05 

.20 

0.46 

1.37 

0.06 


O 
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Table V (continued) 


Phase 

R UMa 

.25 

0.48 

.30 

0.51 

.35 

0.57 

.40 

0.69 

. 45 

0.81 

+0.50 

0.95 


R Vir 

SS Vir 

1.37 

0,10 

1.38 

0.17 

1.42 

- 

1.48 

- 

1.57 

- 

1.65 

- 
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and phase. 

After these eyeball curves were drawn, the 
average 3 . 5y magnitude as a function of phase was 
recorded in table V for phases —0.5 to +0.5 in steps 
of 0.05. It must be emphasized that table V is a mean 
curve; individual cycles may deviate 10 to 20% from the 
magnitude but, hopefully, will retain the general shape 
of the curve. 

There are some general remarks to be made con- 
cerning the 3.5y light curves as plotted in figures 4-7. 
One comment is that there are cycle to cycle variations- 
the 3.5y magnitudes are not identical at similar phase 
points in differing cycles. Part of this alteration 
may be masked within the observational errors , and 
part may be due to slightly changed periods among 
different cycles causing errors in the phase calcul- 
ations. However, where the data are most complete, 
as in x Cyg for example , the cycle to cycle changes 
appear to be real and can vary as much as 10-20% 
from the adopted mean curve. Inspection of the figures 
does not allow any statement about the relation of 
cycle to cycle variations to period or to spectral type. 

A second feature evidenced by the figures concerns 
the shape of the 3.5u light curve. Again, for stars 
with the most complete data, the general shape of the 
3 . 5y curve may be retained, but its magnitude is 
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altered. The shape of the' light curve is not a 
simple trigonometric function and certainly is not 
a reproduction of the visual curve reduced in amplitude 

One of the most significant aspects of the figures 
is the phase lag between the 3 . 5y and visual light 
curves . This lag is the rule rather than the exception 
Maximum brightness at 3 . 5y usually occurs between 
phase 0.1 and 0.2, whereas minimum brightness occurs 
(about one-half cycle later) near phase -0.4 to -0.3. 
There are some exceptions to this minimum phase, 
however, notably T Cep and T Cas which will be 
discussed later. 

Probably the second most interesting aspect of 
the 3 . 5y light curves is the existence of humps in 
the rising branch near phase -0.3 to -0.2. The most 
obvious examples are R Agl , T Cas, U Her,x Cyg and 
U Cyg. This hump may be on the order of the esti- 
mated errors in measurement but appears to be real 
as it is evident over more than one cycle (e.g.y Cyg) . 
This feature then tends to move the primary minimum 
3.5y brightness to earlier phases near -0.4 and 
may then produce a secondary maximum near phase -0.25 
and a secondary minimum near phase -0.15. 

Unfortunately, for some of the stars much of 
the data is not complete enough to ascertain the 
presence of humps. For the objects with good data 
coverage, these secondary extrema do not seem to occur 
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in all cases. Some light curves merely possess an 
inflection point or 1 shape alterations on the rising 
portion of the 3 . 5y light curve near phase -0.25. 
True humps seem to be the rarity while inflections 
are more common. 

Each Mira type variable star seems to have a 
personality all its own and its behavior cannot be 
sufficiently described by a few parameters such as 
mean spectral type at maximum, mean period and mean 
amplitude at distinct wavelengths. Comments per- 
taining to some of the individual stars in 
figures 4, 5 and 6 are included in the following 
section to provide additional information. 

R Boo {M 4 . 5e , 224 d , IRC + 30260}: Cycle 2 
which had its visual maximum near J. D. 2440730 
was brighter than normal and appears to be re- 
flected in the 3.5y data also. 

R Tri {M 4e , 266 d , IRC + 30044}: This star 
demonstrates reasonably good reproducibility at 
visual and 3.5y wavelengths. It was measured 
radiometrically by Pettit and Nicholson (1933) . 

R Vir (M 4 . 5e , 145 d , IRC + 10256}. 

R Aql (M 6.5e, 306 d , IRC + 10406}: The period 

j 

of this star has been decreasing from near 350 , at 
the turn of the century, to 293 d since 1954 (Ahnert, 
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1964, 1965; Schnellar, 1965). Campbell (1955) 
cites a 30 6 d period, whereas a 2 83 d period was 
used in the phase calculations. This object has a 
hump on the ascending portion of the 3.5y light 
curve. It is also a known OH (Wilson and Barrett, 
1970) and H 2 0 microwave emitter (Schwartz and 
Barrett, 1970) . 

o Cet {M 5.5e, 331 d , IRC + 00030}: This 
object had an unusually bright visual maximum near 
J. D. 2440460, or 3 . 5y cycle 2, which also may be 
reflected in the 3.5y data. o Cet is a known binary, 
but its companion should not affect the 3.5y observa- 
tions at all. Unfortunately, the 3 . 5y data quality 
and coverage could be improved for this prototype. 

It was also monitored by Pettit and Nicholson (1933) . 

R Leo {M 7e, 313 d , IRC + 10215} : Another Pettit 
and Nicholson (1933) object. 

X Oph {M 6 . 5e , 334 d , IRC + 10366}: Another 
Pettit and Nicholson (1933) source. It is another 
known binary but the 3 . 5y measurements should not be 
influenced by the K type companion. 

R UMa { M 4 . 5e , 301 d , IRC + 70099} . 

j 

R Aqr {M 7e + pec, 386 , IRC - 20642}: This star, 
with its peculiar spectrum, is extremely low in the 
Minnesota sky such that extinction uncertainties may 
seriously affect the data quality (which is limited 
to near zero phase) . Pettit and Nicholson (1933) 
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also measured this star radiometrically . 

R Cnc {M 6.5e, 362 , IRC + 10185}: Another 
Pettit and Nicholson (1933) source. 

T Cep (M 6.5e, 393 d , IRC + 70168}: This star 

is the exception to the rule at 3.5y. It possesses 
a very small amplitude, but even with the observa- 
tional errors, it seems unique in its phase with 
respect to the visual. This source seems to minimize 
near zero phase and maximize near phase 0.50, in 
direct contrast to the other variables. Lockwood 
and Wing (1971) support this small amplitude and the 
unusual shape at 1.04y. They also comment that the 
period may be decreasing. Campbell (1955) quotes a 
mean period of 39 3, while Keenan (1966) mentions a 
reduction of the period since 1951 to 37575. A 
382 period seemed to fit the recent visual 
behavior and was used in the phase calculations. 

This object has been classified as an S star by 
Spinrad and Newburn (1965) . 

S CrB (M 6 . 5e , 360 d , IRC + 30272}: The 3.5y 

observations repeat themselves quite well. This 
object has been reported as an I^O microwave emitter 
(Schwartz and Barrett, 1970) . 

R LMi {M 7e , 372 d , IRC + 30215}: The 3.5]i data 
coverage is limited for this star also monitored by 
Pettit and Nicholson (1933) . 
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U Ori {M 6 . 5e , 373 d , IRC + 20217): This is an 

OH (Wilson and Barrett, 1970) and an H^O source 
(Schwartz and Barrett, 1970). 

R Ser {M 6 . 5e , 357 d , IRC + 20285}: Classified 
as a possible S star by Spinrad and Newburn (1965) . 

R Aur (M 6.5e, 462 d , IRC + 50141}. 

R Cas {M 7 . 5e , 430 d , IRC + 50484}: This object 
seems to show true cycle to cycle alterations in 
agreement with Lockwood and Wing (1971) . The visual 
maximum near J. D. 2441020 was brighter than normal 
and corresponds to 3 . 5y cycle 3, which also supports 
this fact. 

_3 

T Cas (M 7 . 5e , 445 , IRC + 60009}: This source 
has one of the most pronounced humps on the rising 
branch of the 3.5u light curve. This extreme hump 
is also seen in the mean visual curve (Campbell, 

1955) and in the 1.04y light curve (Lockwood and 
Wing, 1971). This object has also been classified 
as an S type star (Spinrad and Newburn, 1965). 

U Her {M 6.5e, 405 d , IRC + 20298}: This H 2 0 
microwave emitter (Schwartz and Barrett, 1970) seems 
to have an inflection point on the rising portion 
of the 3.5 light curve. 

R And { S4 , 7e , 409 d , IRC + 40009}: Cycle 3, 
which had its visual maximum near J. D. 2441050, 
was visually brighter than normal but this brightness 
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is not evidenced by the 3.5y measurements. 

W And { S8 , 2e , 397 d , IRC + 40037}. 

R Cyg {S3, 9e , 426 d , IRC + 50301}: Another object 
measured by Pettit and Nicholson (1933). The 3.5y 
data reproduction is much better than could be 
expected for the magnitude range of this source. 

It has one of the largest Ul™2) 3.5y variations 
of any of the stars of this program. 

X Cyg { S7 , 2e , 409 d , IRC + 30395}: This star 
has the most complete and continuous 3 . 5y data 
coverage. The visual maximum near J.D. 2440480 was 
a little brighter than average, but does not seem so 
at 3 . 5y in cycle 2. Cycle 3 at 3 . 5y appears dimmer 
than average and may be supported by the visual 
observations. This object has the most well docu- 
mented hump at 3.5y in the ascending branch of the 
light curve. This hump is only slightly reflected 
in the visual curves (Campbell, 1955) but seems 

O 

present at 8500 A (Hetzler, 1936) . This anomaly 
may or may not be present in every cycle at 1.0 4y 
(Lockwood and Wing, 1971). This star has the 
largest visual range (~8 I I l 2) of any well observed 
Mira variable. In the infrared, the 3.5y range is 
definitely not the largest. 

S Cep {C7 4 e, 485 d , IRC + 80048}. 



71 


V CrB {C 6 2 e, 358 d , IRC + 40273). 

U Cyg (C7 2 e 7 463 d , IRC + 50324): This carbon 
star has an inflection point near phase -0.25. 

RS Cyg (C 8 2 e, 419 d , IRC + 40397): This star 
is unusual because it demonstrated no detectable 
variation at 3.5y. This limited alteration is in 
reasonable agreement with the AAVSO visual data. 

SS Vir {C 6 3 e, 356 d , IRC + 00217). 

In order to describe the 3.5u variations of 
Mira type stars, a limited Fourier analysis method 
was employed. The mean magnitudes from Table V for 
each star were converted to relative intensities and 
then fitted by a seven term Fourier series. 

Equation 4 is the series which was used in the 
curve fitting. 

I- _ ( 6 ) = K + A sin 0 + B cos 0 

3 . 5u 

+ C sin 26 + D cos 20 (4) 

+ E sin 30 + F cos 30 

I 3 is the 3 . 5u intensity expressed in relative 

units; 0= 2ttoc, where a is the visual phase; K-F are 
constants to be determined from the analysis. 

The 7 unknowns were calculated by solving 7 
simultaneous equations for intensity as a function of 
phase where the data were most complete (hopefully, 
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over the entire cycle) . Seven phases were chosen 
from phases 0, .16, .25, .33, .50, -.33, -.25, 

-.16. If the data were limited to about one-half 
cycle, the last two terms of equation 4 were omitted. 

The point must be stressed that this analysis 
was not accomplished with the aid of a computer but 
represents a simple attempt to determine the relative 
contributions of the various terms of the series. As 
seen from figures 4-7, the 3.5y variations are 
typically described by a sine curve which lags the 
visual light curves. Table VI lists the program 
stars and the constants determined from this crude 
method. The amplitude in the second column is the 
ratio of the maximum to minimum observed intensity. 
Figure 8 plots the observed and computed relative 
fluxes as functions of phase for three representative 
stars :x Cyg , T Cas and S Cep. This figure is in- 
cluded to demonstrate what shapes of curves may be 
described by this Simple method. 

It must be strongly emphasized, at this time, 
that the previous curve fitting, Fourier method is 
purely descriptive and does not imply that the over- 
tones of the fundamental frequency are present. Mira 
variables are most surely non-linear oscillators and 
the overtones are not simple multiples of the funda- 
mental frequency. 



Table VI 

Fourier Analysis Parameters 


Star 

Spectral 
Type at 
Maximum 

Amplitude 

K 

A 

B 

C 

D 

E 

F 

R And 

S4,7e 

2.59 

1.62 

.72 

.08 

.10 

-.07 



W And 

S8 ,2e 

1.70 

1.40 

.30 

.195 

.005 

-.025 



R Aql 

M6 ,5e 

1.39 

1.11 

.105 

.09 

.055 

-.06 

-.005 

0 

R Aqr 

M7e+pec 

1.43 

1.05 

.115 

.15 

.08 

.015 



R Aur 

M6. 5e 

1.60 

1.21 

.255 

.135 

.06 

-.06 

.015 

.005 

R Boo 

M4.5e 

1.49 

1.20 

.20 

.12 

.04 

0 

.015 

0 

R Cnc 

M6 . 5e 

1.36 

1.08 

.25 

.19 

-.04 

-.005 



R Cas 

M7e 

1.61 

1.23 

.22 

.155 

.045 

-.025 

.02 

-.025 

T Cas 

M7 . 5e 

1.41 

1.18 

.14 

.025 

.08 

-.06 

0 

-.005 

S Cep 

C7 4 e 

1.81 

1.28 

.06 

.36 

.11 

.085 



T Cep 

M6 . 5e 

1.17 

1.07 

.05 

-.06 

-.015 

.01 

0 

-.01 

o Cet 

M5 . 5e 

1.57 

1.32 

. 22 

. 15 

.03 

-.06 

-.01 

.02 

S CrB 

M6.5e 

1.57 

1.26 

.24 

.11 

.01 

-.07 

.01 

-.01 

V CrB 

C6 2 e 

1.99 

1.42 

.25 

.48 

-.01 

.14 

-.04 

.04 

R Cyg 

S3 , 9e 

3.04 

2.16 

1.04 

.42 

.12 

.04 

.015 

.045 


-4 

u» 



Table VI (continued) 


Star 

Spectral 
Type at 
Maximum 

Amplitude 

K 

A 

B 

C 

D 

E 

F 

U Cyg 

C7 2 e 

2.82 

1.78 

.45 

.64 

.15 

-.035 

.08 

.10 

X Cyg 

S7,2e 

1.89 

1.27 

.34 

.21 

.08 

-.04 

.05 

.015 

U Her 

M6 . 5e 

1.50 

1.26 

.12 

.23 

.05 

-.02 

.03 

.02 

R Leo 

M7e 

1.47 

1.24 

.215 

.05 

-.03 

-.07 

.01 

.03 

R LMi 

M7e 

1.36 

1.08 

.155 

.10 

.04 

-.055 



X Oph 

M6.5e 

1.37 

1 « 4 . 6 . 

.10 

.115 

.02 

-.05 

-.01 

-.005 

U Ori 

M6.5e 

2.06 

1.38 

.415 

.215 

.11 

-.045 

.015 

-.045 

R Ser 

M6. 5e 

1.40 

1.18 

.15 

.12 

0 

-.045 

-.01 

.015 

R Tri 

M4e 

1.33 

1.07 

.14 

.035 

.03 

-.025 

.01 

-.01 

R UMa 

M4 . 5e 

1.64 

1.34 

.22 

.20 

-.035 

-.06 

.005 

.03 

R Vir 

M4 . 5e 

1.31 

.98 

.28 

.04 

-.005 

.01 
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PHASE 

Figure 8. A comparison of the mean 3.5u relative flux 
versus visual phase and the calculated flux 
from the limited Fourier curve fitting 
analysis . 
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VI 

OTHER WAVELENGTH OBSERVATIONS 

In addition to the 3.5y observations, measure- 
ments at other wavelengths were obtained on a limited 
number of the same program stars. These data were 
collected with the Minnesota photometric systems used 
on the O’Brien telescope in Minnesota and on the 
Mt. Lemmon infrared telescope in Arizona. The details 
of the photometer and the absolute calibration for 
the various wavelengths have been previously described. 
The method of data acquisition and analysis was 
similar to that prescribed for the 3.5y observations 
in preceding sections. 

Table VII lists some of the calibration stars 
by their HR number, name, spectral type, visual 
magnitude, and magnitudes at other effective wave- 
lengths. Most of the listed magnitudes are quoted 
to an accuracy + 10%. 

Appendix 3 tabulates the observational data for 
the program stars at wavelengths in addition to 3.5y. 

The results are again quoted to an accuracy of + 10%. 

The data are by no means complete over entire cycles 
because another observer (Forrest, 1973) is monitoring 
many of these same program stars, primarily at effective 



Table VII 

Magnitudes of Calibration Stars 
Spectral* Magnitudes at Effective Wavelength (microns) 


HR 

Name 

Type 

m * 
v 

2.3 

3.5 

4.8 

8.6 

10.7 

12.2 

18 

337 

6 And 

MO III 

2.03 

-1.8 

-2.05 

-1.7 

-2.0 

-2.1 

-2.1 

- 

1457 

a Tau 

K5 III 

0.86 

'-2.9 

-3.02 

-2.9 

-3.0 

-3.1 

-3.1 

-3.1 

2061 

a Ori 

M2 lab 

0.80 

-4.20 

-4.54 

-4.20 

-4.75 

-5.55 

-5.50 

-5.65 

2990 

6 Gem 

KO III 

1.15 

-1.1 

-1.23 

-1.1 

-1.2 

-1.2 

-1.2 

-1.3 

4069 

y UMa 

MO III 

3.04 

-0.8 

-1.00 

-0.8 

-1.0 

-1.1 

-1.1 

- 

5340 

a Boo 

K2 IIIp 

0.06 

-3.02 

-3.12 

-2.95 

-3.16 

-3.20 

-3.20 

-3.2 

6406-7 

a Her 

M5 II 

3.1 

-3.44 

-3.70 

-3.48 

-3.79 

-4.05 

-4.00 

-4.1 

7001 

a Lyra 

AO V 

0.04 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

8316 

U Cep 

M2 la 

3.99 

-1.8 

-2.2 3 

-2.2 

-3.1 

-4.3 

-4.1 

-4.7 

8775 

6 Peg 

M2 II-III 

2.56 

-2.24 

-2.45 

-2.27 

-2.45 

-2.55 

-2.52 

-2.6 


* Hoffleit (1964) . 


--j 
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wavelengths of 3.5, 4.9, 8.4 and 11.0 microns. 

Many of the stars of this research project have 
been measured at other wavelengths of 3.5, 4.9, 8.4, 
and 11. Oy by Gillett, Merrill and Stein (1971), 

Gehrz and Woolf (1971), and Dyck, Forrest, Gillett, 

Stein, Gehrz , Woolf and Shawl (1971). These measure- 
ments have been summarized in table VIII in order to 
provide a more complete description of the energy 
distribution of long period variables. It must be 
stressed that the magnitudes quoted in table VIII 
have no phases listed. The date of the observation 
is not available nor published in the articles 
cited. Therefore, also included in columns 6 and 
7 of the table are the observed magnitude ranges at 
3 . 5y as taken from table IV and appendix I. Column 8 
which lists the observed [3.5y]-[lly] color index is 
included to provide an estimate of the excess radiation 
at lip arising from the circumstellar shells. 

A multitude of authors have published observational 
data at wavelengths shorter than 3 . 5y for many of these 
same stars (Mendoza and Johnson, 1964; Johnson, Mitchell, 
Iriarte and Wisniewski, 1966). From these shorter 
wavelength magnitudes and, in conjunction with the 
longer wavelength data, typical spectral energy 
distributions may be constructed for some of the 
representative LPV stars. Two examples are the carbon 



Table VIII 

Summary of 3.5, 4.9, 8.4 and 11 . Op Magnitudes of the Program Stars 

Magnitude at effective wavelength (microns) 

3 . 5p 


Star 

3.5 

4.9 

8.4 

11.0 

max 

min 

[3. 5y ] — [lly] 

R Aql 

-1.19 

-1.24 

-1.76 

-2.87 + 

-1.5 

-1.1 

1.7 

R Boo 

1.18 



0.42 + 

1.1 

1.7 

0.8 


1.26 

1.10 

0.71 

0.26 : 



1.0 


1.34 

1.03 

0.64 

0.10 * 



1.2 

R Cnc 

-1.33 

-0.88 

-1.49 

-2.55 + 

-1.3 

-1.0 

1.2 


-1.16 

-1.51 

-1.94 

-2.56 * 



1.4 

T Cnc 

0.42 

0.59 

-0.56 

-0.65 * 

0.2 

0.4 

1.1 

R Cas 

-2.06 

-2.18 


-4.08 + 

-2.6 

-2.1 

2.0 


-2.00 

-2.43 

-3.06 

-4.10 * 



2.1 

T Cas 

-1.52 



-2.61 + 

-1.8 

-1.4 

1.1 


-1.58 

-1.67 

-2.22 

-2.94 * 



1.4 

WZ Cas 

0.46 

0.54 

0.23 

-0.04 * 

-0.1 

0.4 

0.5 

S Cep 

-1.64 

-2.02 


-2.91 + 

-1.8 

-1.1 

1.3 


-1.47 

-1.73 

-2.63 

-3.11 * 



1.6 

T Cep 

-2.22 

-2.12 

-2.72 

-3.15 * 

-2.4 

-2.2 

0.9 

o Cet 

-3.49 



-5.45 + 

-3.7 

-2.9 

2.0 


-3.49 

-3.81 

-4.59 

-5.63 : 



2.1 

V CrB 

1.29 

0.69 

-0.11 

-0.85 * 

0.3 

1.3 

2.1 


-j 

VD 



Table VIII (continued) 


Magnitude at effective wavelength (microns) 

3 . 5y 


Star 

3.5 

4.9 

8.4 

11.0 

max 

min 

[3 . 5uJ - 

S CrB 

-0.60 

-0.98 

-1.98 

-2.83 + 

-0.8 

-0.3 

2.2 


-0.74 

-1.19 

-2.18 

-3.12 * 



2.4 

R Cyg 

0.20 

0.19 

-0.68 

-1.27 * 

0.2 

1.4 

1.5 

U Cyg 

0.65 

0.24 

-0.57 

-1.40 + 

-0.1 

1.0 

2.1 


0.00 

-0.28 

-1.00 

-1.61 * 



1.6 

RS Cyg 

1.19 

1.31 

0.66 

0.52 * 

0.8 

1.1 

0.7 

X Cyg 

-2.25 

-2.57 

-3.21 

-4.00 + 

-3.2 

-2.4 

1.7 


-2.95 

-3.06 

-3.35 

-4.16 * 



1.2 

RY Dra 

-0.37 

0.03 

-1.10 

-1.20 * 

-0.4 

-0.4 

0.8 

U Her 

-0.78 

-1.11 

-1.67 

-2.59 + 

-1.0 

-0.6 

1.8 

R Leo 

-3.25 

-3.59 

-4.02 

-4.93 + 

-3.4 

-2.9 

1.7 


-3.17 

-3.39 

-3.92 

-4.65 : 



1.5 


-3.21 

-3.41 

-3.90 

-4.56 * 



1.3 

R LMi 

-1.15 

-0.96 

-2.15 

-2.83 + 

-1.4 

-1.0 

1.7 


-1.14 

-1.60 

-2.11 

-2.82 * 



1.7 

X Oph 

-1.45 

-1.49 


-2.76 + 

-1.6 

-1.0 

1.3 


00 

o 



Table VIII (continued) 


Magnitude at effective wavelength (microns) 

3.5 y 


Star 

3.5 

4.9 

8.4 

11.0 


max 

min 

[3.5y]- [lly] 

U Ori 

-1.20 

-0.93 

-1.80 

-2.82 

+ 

-1.7 

-1.0 

1.6 


-1.37 

-1.66 

-2.11 

-3.00 

* 



1.6 

W Ori 

-0.72 

-0.33 

-1.24 

-1.74 

* 

-1.1 

-0.8 

1.0 

R Ser 

0.18 

-0.03 

-0.58 

-1.26 

* 

-0.4 

0.0 

1.4 

R Vir 

1.49 

1.08 

0.78 

0.64 

* 

1.3 

1.6 

0.8 

SS Vir 

0.17 

0.46 

-0.56 

-0.99 

* 

0.0 

0.3 

1.2 

: Dyck et 

al. (1971). 







+ Gehrz and Woolf 

(1971) . 







* Gillett 

, Merrill 

and Stein 

(1971) . 
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star U Cyg near infrared maximum and the mild S star 

X Cyg near infrared (3.5p) minimum, as plotted in 

2 

figure 9 . The figures depict AF^ in watts/cm versus 
wavelength. U Cyg has been brightened by 4 magnitudes 
in the graph. The short wavelength U - L points 
are from Mendoza and Johnson (1965) . The longer 
wavelengths are from Gehrz and Woolf (1971) , and 
Gillett, Merrill and Stein (1971). These two 
figures are included to emphasize the differences 
in spectral shape among the carbon rich and the 
oxygen rich stars. These differences will also be 
discussed in a later section. 

In addition to these multi-wavelength measure- 
ments, brightness changes of long period variable 
stars at single selected wavelengths have been 
monitored by individuals and groups. In the visual 
region, the American Association of Variable Star 
Observers (AAVSO) monitors many LPV stars. Near 

O 

8500 A, Hetzler (1936) catalogued the variations. 
Landolt measures UBV of many Mira variables . Pettit 
and Nicholson (1933) were the first to measure the 
radiometric changes. More recently, Lockwood and 
Wing (1971) and Lockwood (1972) have measured changes 
at wavelengths near one micron. 

Lockwood and Wing (1971) produced narrow band 
photometry on a group of M type Miras at wavelengths 
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BVRIJ K LM 0 



WAVELENGTH (MICRONS) 

Figure 9. Relative energy distribution (AF^ in 
watts/cm^ versus wavelength) for' the 
mild S star (x Cyg) and a carbon star 
(U Cyg) . U Cyg has been brightened by 
4 magnitudes. 


C/ 
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of .78, .87, .88, 1.04 and 1.05 microns. In the 
oxygen rich stars, they believe that the 1.04y 
magnitudes are uncontaminated by absorption, and, 
therefore, represent the true continuum point of the 
photospheric temperature. They have measured light 
curves in 5 colors for two stars (o Cet and T Cas) . 

For many other variables they produced light curves 
at 1.0 4u only. 

By pure chance, but also because of apparent 
brightness, Lockwood and Wing (1971) measured many 
of the same M stars which comprise this project (or 
vice versa) . Table IX lists the LPV stars common 
on the two projects, their spectral type as determined 
by Lockwood and Wing (1971) and Lockwood (1972) , the 
mean visual range (Campbell, 1955), the range in 
magnitude at 8500 A (Hetzler, 1936), and the 1.04 m 
and 3.5u magnitude variations. It merits attention 
that the observed extrema may not be indicative of 
the true range of stellar variability because of the 
limited 3.5y data coverage. 

The 1.04y data of Lockwood and Wing (1971) and 
Lockwood (1972) were graphed as a function of visual 
phase in the same manner as was done for the 3.5p 
observations of this project. A curve was then drawn 
through the individual points "by eye". The shape of 
the curve was determined from the magnitudes on 
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Table IX 

Range of Mean Magnitudes at Selected Wavelengths 


Star 

Spectral 

Type 


Visual* 

8500A+ 

1.04ii 3 . 5y 

R And 

S4 , 7e 

a 

7™4 

3™6 

- 

1™2 

W And 

S8 ,2e 

a 

6.3 

- 

- 

0.6 

R Aql 

M6.0-M8.5 

m 

• 

5.4 

2.6 

0.9 

0.4 

R Aur 

M6.0-M9.5 

• 

5.8 

- 

1.2 

l 0 * 5 

R Boo 

M2-M8.2 

• 

5.1 

- 

0.75 

0.7 

R Cnc 

M6 . 5e 

a 

4.4 

- 

- 

- 

R Cas 

M6 . 0-M10 

m 

9 

5.6 

2.9 

1.65 

>_0.5 

T Cas 

M6.0-M9.0 

9 

9 

4.0 

- 

0.7 

>_0 . 4 

T Cep 

M5.5-M8. 8 

m 

9 

4.3 

2.2 

0.55 

0.2 

o Cet 

M4 . 5-M9 . 2 

* 

5.7 

2.8 

1.1 

0.8 

S CrB 

M6 . 5e 

a 

5.6 

- 

- 

0.5 

R Cyg 

S3 ,9e 

a 

6.4 

- 

1.55 

1 . 2 

X Cyg 

M2-M9.0 


8.2 

4.0 

2.0 

0.7 

U Her 

M7 . 5-M9.5 


5.0 

- 

1.35 

>0.5 

R Leo 

M6.0-M9.2 


4.2 

2.4 

1.1 

>0.5 

R LMi 

M6 . 5-M9 . 0 


5.6 

- 

0.95 

>0.4 

X Oph 

M7.0-M8.5 

m 

9 

2.0 

- 

0.55 

0.4 

U Ori 

M7.5-M9.0 

w 

9 

5.7 

- 

1.2 

0.7 

R Ser 

M6 . 5e 

a 

6.5 

- 

- 

>0.4 

R Tri 

M3-M8 . 5 

• 

» 

5.5 

2.5 

0.85: 

0.5 

R Vir 

M4 . 5e 

a 

4.6 

1.6 

- 

0.2 

Mean 



5.4 

2.7 

1.1 

0.7 


* Campbell (1955) . 

+ Hetzler (1936) . 

: Lockwood and Wing (1971) . 
a Keenan (1966 ) . 
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individual cycles; the mean value was deduced from the 
constraints imposed by the different cycles. From 
these averaged curves, a mean 1.04y magnitude as a 
function of phase was determined and is listed in 
table X. Figure 10 is a plot of the mean magnitudes 
as functions of phase for the visual, 1.04y and 3.5y 
wavelengths of the representative stars x C Y9» R Cyg , 
and T Cas . This figure is included to depict the 
relative shapes and amount of variation at the different 
wavelengths. The mean visual curves are from Campbell 
(1955) . 

Because the 1.0 4y and 3.5y magnitudes are 
supposedly free of absorptions and represent a true 
continuum level of the stellar radiation, a color 
index of these two magnitudes should provide a good 
temperature indicator of the photosphere. Armed 
with this color index as a function of phase, the 
parameters of temperature, angular diameter, luminosity 
and assorted others seem to be the next quantities 
which should be computed. 



Table X 

Adopted 1.04y Magnitudes as a Function of Phase 


Phase 

R And 

W And 

R Aql 

R Aur 

R Boo 

R Cnc 

R Cas 

0.00 

2.95 

- 

0.81 

1.25 

3.54 

1.06 

0.00 

.05 

2.48 

- 

.80 

1.22 

3.44 

1.06 

-0.02 

.10 

2.00 

- 

.85 

1.23 

3.36 

1.10 

0.02 

.15 

- 

- 

.95 

1.29 

3.35 

1.20 

0.12 

.20 

- 

- 

1.07 

1.41 

3.39 

1.30 

.31 

.25 

- 

- 

1.21 

1.61 

3.47 

1.39 

.50 

.30 

- 

2.40 

1.35 

1.79 

3.56 

- 

.68 

.35 

- 

2.62 

1.51 

2.02 

3.67 

- 

.85 

.40 

- 

2.84 

1.63 

2.31 

3.79 

- 

1.04 

+0.45 

- 

3.09 

1.70 

2.37 

3.96 

- 

1.26 

+0.50 

- 

3.26 

1.74 

2.36 

4.03 

- 

1.50 

-0.45 

4.40 

3.38 

1.66 

2.19 

4.12 

- 

1.69 

.40 

4.72 

3.39 

1.50 

1.97 

4.15 

- 

1.71 

.35 

4.72 

3.20 

1.35 

1.79 

4.07 

- 

1.56 

. 30 

4.64 

2.69 

1.23 

1.67 

3.96 

- 

1.32 

. 25 

4.54 

— 

1.17 

1.57 

3.87 

- 

1.04 


CO 



Table X (continued) 


Phase 

R And 

W And 

R Aql 

R Aur 

R Boo 

R Cnc 

R Cas 

.20 

4.40 

- 

1.15 

1.52 

3.80 

+1.34 

0.74 

.15 

4.24 

- 

1.12 

1.46 

3.75 

1.29 

0.42 

.10 

3.94 

- 

1.03 

1.39 

3.70 

1.21 

0.18 

-0.05 

3.43 

- 

0.91 

1.32 

3.63 

1.12 

0.04 

0.00 

2.95 

- 

0.81 

1.25 

3.54 

1.06 

0.00 


Phase 

T Cas 

T Cep 

o Cet 

S CrB 

R Cyg 

X Cyg 

0.00 

0.85 

0.09 

-1.26 

- 

3.03 

-0.02 

.05 

.81 

.11 

-1.34 

1.57 

2.97 

-0.06 

.10 

.86 

.16 

-1.36 

1.64 

2.97 

-0.03 

.15 

.98 

.22 

-1.34 

1.72 

3.02 

0.08 

.20 

1.10 

.28 

-1.28 

1.83 

3.14 

.24 

.25 

1.24 

. 31 

-1.21 

1.94 

3.28 

.46 

.30 

1.39 

. 33 

-1.12 

2.05 

3.47 

.71 

. 35 

1.53 

. 37 

-1.00 

- 

3.66 

.95 

.40 

1.60 

.44 

-0.86 

- 

3.87 

1.21 

+0.45 

1.58 

.52 

-0.71 

— 

4.08 

1.55 


00 

oo 



Table X (continued) 


Phase 

T Cas 

T Cep 

o Cet 

S CrB 

R Cyg 

X Cyg 

+0.50 

1.50 

.60 

-0.52 

- 

4.27 

1.75 

in 

o 

1 

1.37 

.44 

-0.36 

- 

4.44 

1.77 

. 40 

1.26 

.29 

-0.28 

- 

4.58 

1.61 

.35 

1.14 

.21 

-0.26 

- 

4.66 

1.30 

.30 

1.05 

.20 

-0.29 

- 

4.54 

1.04 

.25 

i— 1 

o 

» 

1— 1 

.22 

-0.39 

- 

4.25 

0.79 

.20 

■ 1.08 

.23 

-0.49 

- 

4.02 

0.59 

.15 

1.14 

.22 

-0.57 

- 

3.78 

0.44 

.10 

1.12 

.17 

-0.84 

- 

3.54 

0.31 

-0.05 

1.01 

.11 

-1.10 

- 

3.28 

0.10 

0.00 

0.85 

0.09 

-1.26 

— 

3.03 

-0.02 


Phase 

U Her 

R Leo 

R LMi 

X Oph 

U Ori 

H Ser 

R Tri 

R Vir 

0.00 

1.38 

-0.92 

1.18 

0.84 

1.07 

- 

2.40 

3.64 

0.05 

1.32 

-0.91 

1.16 

0.84 

0.99 

- 

2.33 

3.53 

.10 

1.35 

-0.84 

1.19 

0.87 

0.92 

— 

2.30 

3.42 


00 



Table X (continued) 


Phase 

U Her 

R Leo 

R LMi 

X Oph 

U Ori 

R Ser 

R Tri 

R Vir 

.15 

1.50 

-0.73 

1.26 

0.93 

0.92 

- 

2.30 

3.35 

.20 

1.64 

-0.61 

1.37 

1.00 

0.94 

- 

2.33 

3.34 

.25 

1.77 

-0.48 

1.50 

1.08 

1.05 

2.07 

2.41 

3.38 

.30 

1.90 

-0.36 

1.68 

1.16 

1.19 

2.19 

2.54 

3.44 

.35 

2.05 

-0.24 

1.85 

1.22 

1.36 

2.35 

2.68 

3.53 

.40 

2.20 

-0.11 

2.02 

1.29 

1.54 

2.56 

2.82 

3.65 

0.45 

2.33 

0.02 

2.13 

1.35 

1.73 

2.75 

2.96 

3.79 

±0-50 

2.56 

0.13 

2.14 

1.41 

1.90 

2.90 

3.09 

3.88 

-0.45 

2.71 

0.18 

2.11 

1.38 

2.03 

2.96 

3.14 

- 

.40 

2.74 

0.12 

2.04 

1.30 

2.14 

2.92 

3.13 

- 

.35 

2.69 

0.00 

1.96 

1.21 

2.12 

2.76 

3.06 

- 

. 30 

2.57 

-0.14 

1.87 

1.14 

1.94 

- 

2.97 

- 

.25 

2.42 

-0.27 

1.74 

1.08 

1.80 

- 

2.90 

- 

.20 

2.32 

-0.45 

1.60 

1.03 

1.65 

- 

2.83 

- 

.15 

2.21 

-0.61 

1.47 

1.01 

1.50 

- 

2.73 

- 

.10 

1.99 

-0.73 

1.36 

0.95 

1.34 

- 

2.62 

- 

-0.05 

1.58 

-0.84 

1.26 

0.89 

1.20 

- 

2.50 

- 

0.00 

1.38 

-0.92 

1.18 

0.84 

1.07 

— 

2.40 

3.64 



MEAN 3.5/1 MAG. MEAN 1.04/t MAG. MEAN VISUAL MAG. 
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versus visual phase for the mild S star 
X Cyg. Of interest is the "hump" in the 
3.5y curve which may or may not be present 
at other wavelengths . 




PHASE 

Figure 10b. Same as figure 10a but for the S star 
R Cyg. Note the absence of the hump 
at all wavelengths on the ascending 
portion of the light curve. 





MEAN 3.5 ft MAG. MEAN 1.04 ft MAG. MEAN VISUAL MAG 



PHASE 

Figure 10c. Similar to figure 10a but for the M star 
T Cas . Note the presence of a hump at 
all wavelengths on the ascending branch 
of the light curve. 
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VII 

TEMPERATURES, ANGULAR SIZES AND LUMINOSITIES 

A multitude of attempts have been made to 
calculate the effective or color temperatures of 
the photospheres of M giant stars. One of the 
earliest attacks on this problem was by Pettit and 
Nicholson (1933) who used thermocouples to measure 
the radiometric magnitude and then derive temperatures 
for a few long period variable stars. Smak (1964) 
used spectroscopic and photometric UBV data to 
derive assorted parameters for red giants . Johnson 
(1966) used UBVRIJKLMN photometry to calculate 
effective temperatures for main sequence stars and 
cool giants. Lee (1970) employed much the same 
procedure as Johnson and calculated the effective 
temperatures as functions of spectral type for both 
M giant and supergiant stars. Smak (1964) covered 
the spectral range MO III - M8 III, Johnson (1966) 
did the same for G5III-M6III and Lee did so for MOIII- 
M6 III. 

Near minimum light, some of the stars of this 
present study reach spectral types later than M8 for 
which no effective temperatures are available. An 
extension of the spectral type-temperature scale was 


needed. 
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Two alternative methods of color temperature 
determinations are both based on the observational 
data. Both methods require the formation of a color 
index (C.I.), which is defined as the difference 
between the 1.0 4y and 3 . 5y magnitudes and is given 
by equation 5. 

COLOR INDEX = C.I. = [ 1 . 04y ] - [ 3 . 5y ] (5) 

The calculation of a color temperature from a 
color index involves the major assumption that the 
star radiates like a black body at the two respective 
wavelengths. In other words, the fluxes at each 
wavelength represent radiation emanating from the 
same temperature or optical depth. This determination 
of temperature also assumes that the opacities at each 
wavelength alter in the same manner throughout the 
stellar cycle. 

As a first attempt to calculate temperatures, 
a color index-spectral type relation must be deter- 
mined. This method utilizes the 1.04y observations 
of representative late type giant stars by Lockwood 
(1972) in conjunction with the 3.4y and 3 . 5y data 
on the same stars as gleaned from the literature and 
the present study. The ensuing problem of conversion 
of spectral types into temperatures and the extension 
of the color index to later spectral types are not 
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simple tasks. 

Table XI lists the name, spectral type, 1.04y 
magnitude of Lockwood (1972), the 3 . 4y or 3 . 5y 
magnitudes from various sources, the 3.5p mean magni- 
tude and the [1 . 04y ] - (3 . 5y } color index for late type 
giant stars. 

It deserves to be mentioned at this point that 
Johnson's and Lee's effective wavelengths were 3.4]j 
and not 3 . 5y as in the other observations listed. 
Outside of observational uncertainties, there appears 
to be no systematic difference in magnitude dependent 
on effective wavelength for these standard stars. 

As a result, all available magnitudes were used to 
determine the weighted mean magnitude for the 
individual stars of table XI. 

The next step in the temperature determination 
is to assume or derive a temperature-spectral type 
relation. Table XII lists the effective temperatures 
for red giant spectral types as determined by various 
authors. Also included in the last column, as 
summarized from table XI, are the mean color index 
and, in parentheses, the number of stars used to 
determine this index. 

There are two major problems inherent in the 
color index-effective temperature method as described 
above. First, the non-variable red giants of later 



Table XI 

Observational Summary 


Star 

Spectral 

Type 

1 . 0 4p a 
Mag 

3 . 5y b 
Mag 

3 . 4p c 
Mag 

B And 

MO III 

-0.40 

-2.05 

-2.07 

P UMa 

MO III 

0.60 

-1.00 

- 

75 Leo 

MO Illb 

2.79 

- 

- 

a Vul 

MO III 

2.00 

0.42 

- 

y Eri 

Ml III 

0.53 

- 

-0.97 

X Peg 

M2 III 

2.10 

- 

0.46 

a Cet 

M2 III 

-0.19 

-1.91 

-1.76 

tt Leo 

M2 III 

1.99 

- 

- 

83 UMa 

M2 I I lab 

1.91 

0.25 

- 

S Peg 

M2+IIb-IIa 

-0.64 

-2.45 

-2.43 

55 Peg 

M2 III 

2.05 

0.46 

- 

P Gem 

M3 Ilia 

-0.29 

-2.10 

-2.03 

p Per 

M4 Ilia 

-0.31 

-2.20 

- 

to Vir 

M4-M4.5 III 

1.46 

- 

- 

HR 5299 

M4-M4.5 III 

1.34 

- 

-0.53 


of M Giants 


Color 





Mean 

Index 

3.4 p 

3. 5p 

3.5y 

3.5p 

[ 1 . 04p ] 

Mag 

Mag 

Mag 

Mag 

[3.5p] 

-1.97 

-2.05 

- 

-2^04 

1 .64 

-1.04 

-0.95 

- 

-1.00 

1.60 

- 

1.19 

- 

1.19 

1.60 

- 

- 

- 

0.42 

1.58 

H 7 

• 

O 

O 

- 

- 

-0.98 

1.51 

0.61 

- 

- 

0.53 

1.57 

-1.74 

-1.89 

-1.81 

-1.82 

1.63 

- 

0.21 

- 

0.21 

1.78 

- 

- 

- 

0.25 

1.66 

l— 1 

• 

CN 

1 

-2.45 

-2.22 

-2.44 

1.80 

- 

- 

- 

0.46 

1.59 


-2.01 

_ 

-2.05 

1.76 



-2.14 

-0.33 


-2.15 1.84 
-0.33 1.79 
-0.53 1.87 


\D 

-vl 



Table XI {continued) 


Color 


Star 

Spectral 

Type 

1.04u a 

Mag 

3 . 5y b 
Mag 

3 . 4y C 
Mag 

3 . 4 
Mag 

3 . 5p e 
Mag 

3.5u f 

Mag 

Mean 
3 . 5y 
Mag 

Index 
[1 , 04y] - 
[3. 5y] 

VY Leo 

M5 III 

0.98 

- 

-0.99 

- 

-0.96 

- 

-0.98 

1.96 

RR UMi 

M5 III 

0.73 

- 

-1.12 

- 

- 

-0.93 

-1.03 

1.76 

4 

i Ser 

M5 Ilb-IIIa 

0.98 

- 

- 

- 

-1.25 

- 

-1.25 

2.23 

OP Her 

M5 Ilb-IIIa 

1.87 

- 

-0.15 

- 

0.03 

- 

-0.06 

1.93 

R Lyr 

M.5 III 

-0.32 

-2.28 

- 

- 

-2.32 

-2.23 

-2.28 

1.96 

W Cyg 

M4e-M6 

0.47 

-1.70 

- 

- 

- 

- 

-1.70 

2.17 

RZ Ari 

M6 III 

0.78 

- 

-1.24 

- 

- 

- 

-1.24 

2.02 

RS Cnc 

M6 Ilia 

0.20 

-1.94 

-2.00 

- 

-2.07 

-1.82 

-1.96 

2.16 

X Her 

M6 e 

0.42 

- 

- 

- 

-1.65 

-1.64 

-1.65 

2.07 

g Her 

M6 III 

-0.21 

-2.27 

-2.22 

- 

-2.18 

-2.40 

-2.28 

2.07 

CH Cyg 

M6 

1.20 

-1.31 

- 

- 

- 

-1.25 

-1.28 

2.48 

SW Vir 

M7 III 

0.23 

- 

- 

- 

- 

-2.00 

-2.00 

2.23 

RX Boo 

M8 : 

0.35 

-2 . 30 



-2.25 

-2.33 

-2.29 

2.64 


e Gehrz & Woolf (1971). 


a Lockwood (1972). 
b This study. 


c Lee (19 70) . 
d Johnson (1964) . 


f Gillett ,Merrill & 
Stein (1971) . 


00 



Table XII 

Comparison of Effective Temperatures 
and Color Indices as Functions of Spectral Type 


Spectral Type Effective Temperature °K 




Johnson 

(1966) 

Lee 

(1970) 

Smak 

(1964) 

Mean 

Color Index 

MO 

III 

3660 

3700 

3340 

1™61 (4) 

Ml 

III 

3600 

3650 

3200 

1.51 (1) 

M2 

III 

3500 

3550 

3090 

1.67 (6) 

M3 

III 

3300 

3350 

2980 

1.76 (1) 

M4 

III 

3100 

3150 

2850 

1.83(3) 

MS 

III 

2950 

2950 

2710 

1.97(5) 

M6 

III 

2800 

2750 

2500 

2.08 (4) 

M7 

III 



2250 

2.23(1) 

M8 

III 



2000 

2.64 (1) 
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spectral types (M7 and M8) are scarce and ones of 
even later type (M9 and M10) are totally absent in 
the data. Therefore, an extension of the spectral 
type-color index relation to the M9 and M1Q spectral 
types needed for the Miras near minimum cannot be 
directly made. Secondly, the color index method 
directly determines a color temperature and not an 
effective temperature. The conversion of color 

t 

temperature into effective temperature involves 
many assumptions implicit in the calculations. An 
effective temperature is defined as the temperature 
obtained on the hypothesis that the total energy 
output of the star may be represented by one 
temperature. In other words, the luminosity, size 
and distance of the star must be known. A color 
temperature rests on the assumption that the stellar 
photosphere radiates like a black body of one 
temperature at two or more wavelengths. The 
color and effective temperatures may not be coincidental 
because of emission, absorption, limb darkening and 
wavelength dependent opacity. 

The alternative path to the determination of a 
temperature scale avoids the difficulties previously 
mentioned and directly computes a color temperature 
from the measured color index. The 1.04u measurements 
of Lockwood and Wing (1971) and Lockwood (1972) , as 
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well as the 3 . 5u observations of this program, are 
quoted with respect to Vega. This AOV star is defined 
to be zero magnitude at both wavelengths and its 
resulting color index is therefore zero. The 
remaining obstacle was to choose a color temperature 
for this star, a Lyra was defined to have a 10,000, °K 
color temperature. 

Table XIII lists the black body color index in 
magnitudes as a function of temperature. When the 
temperature is on the order of one-third the reference 
temperature or less, the color index varies linearly 
with 1/T. Figure 11 graphs the color index in 
magnitudes as a function of 1/T for the temperatures 
and color indices of red giant stars. For these 
stars, the color index is greater than about 2?0 . 

In a color index range of 2 to 5 , the linear approxi- 
mation between color index and 1/T can be used to an 
accuracy better than about 1%. 

Equation 6 describes the linear fit in the color 
index and 1/T. The constants a and b are determined 
from the graph. 

C.I. = a (1/T c ) + b (6) 

For color indices between 2 and 5 magnitudes , the 
linear approximation requires a=9850 °K and 
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Table XI ii 

[1.04y]- [3.5y] Color Index 
as a Function of Temperature 


Temperature 

Color Index 
Magnitudes 

20,000 °K 

-0™32 

15,000 

-0.22 

10,000 

0.00 

9,000 

0.08 

8,000 

0.18 

7,000 

0.30 

6,000 

0.48 

5,000 

0.74 

4,000 

1.16 

3,000 

1.91 

2,800 

2.12 

2,600 

2.38 

2,400 

2.68 

2,200 

3.06 

2,000 

3.50 

1,900 

3.76 

1,800 

4.05 

1,700 

4.39 

1,600 

4.76 

1,500 

5.18 



COLOR INDEX [l.04p]-[3.5p] 
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Figure 11. [1 . 04u] ~ [ 3 . 5u ] color index as a function 

of (color temperature) - ^. 
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b= -1.41. Solution of equation 6 for the temperature 
results in equation 7 for 5 m > C.I. > 2 m . 


9850 °K 

c C.I. + 1.41 


(7) 


Due to the difficulties of transformation of 
color to effective temperature, and because the color 
indices for the latest spectral types of variable 
red giants are not available, the decision was made 
to employ this second method of color temperature 
determination . 

One point must be stressed: interstellar 
extinction will influence the color temperatures 
deduced from the color index. If this reddening 
follows an approximate A 1 law, the 1.0 4y magnitudes 
will be affected only about half as much as the 
visual magnitudes, while the 3.5]i extinction will 
be less than 20% of the visual extinction. The 
intrinsic colors of the Miras are not well determined 
and corrections to the observed colors stand on an 
unfirm footing. On the basis of available data, 
interstellar extinction has been totally ignored. 

Fortunately, most of the Mira variables in this 
project are the brighter ones, both visually and 
at 3 . 5y , and, therefore, are likely to be nearby 
stars. Also, many are at intermediate or high 
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galactic latitudes and probably suffer little 
extinction. As for the standard stars of table XI, 
they are, again, the brighter ones. The scatter 
of the color indices for these standards, to some 
extent, may reflect the ignorance of the inter- 
stellar extinction. 

Other factors which may influence the color 
temperature determinations are the unknown weak 
bands which affect the 1.04]i magnitudes of S 
stars (Lockwood and Wing, 1971; Lockwood, 1972). 

For this reason, the S stars possess more unreliable 
temperatures and they must be referred to as 
equivalent temperatures. Unfortunately, temperatures 
for carbon stars may not be inferred. The present 
data lack a second wavelength representing a true 
continuum point such that color-indices may be 
calculated . 

Armed with the assumed color index versus 
temperature relation as determined by equation 7 , 
the color temperatures as function of phase through- 
out the stellar cycle may be computed. Once the 
temperature has been determined, other quantities 
such as angular diameter, apparent bolometric magnitude, 
apparent visual magnitude, and titanium oxide 
blanketing may be calculated. 

As a prelude to the actual equations, a brief 
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outline of the procedure employed to determine the 
temperature, angular size, luminosity and visual 
magnitude will be listed. Given the 1.0 4y and 3.5y 
measurements , a color index is formed and thence 
the temperature. The 3 . 5y magnitude and absolute 
calibration give the measured flux, whereas the 
temperature and Planck function result in the 
intensity at 3.5y. From these two quantities the 
solid angle of the source may be deduced. The 
temperature and solid angle are then utilized to 
form the bolometric and visual magnitudes. 

The color index (C.I.) (equation 5) as a 
function of phase was computed from the mean 1.04y 
and 3.5y curves of tables X and V, respectively. 

The color temperatures were then calculated via 
equation 7. The observed flux (F_ c (A)) at a 
given wavelength is determined from equation 8, 

F OBS U) = F 0 U) 10 ~°' 4 [m OBS U) " m 0 (X)1 (8) 

where F Q (A) is the zero magnitude flux for a specific 
wavelength, m Qgg(A) is the observed magnitude at the 
same wavelength, and m o (A) is the zero magnitude 
reference. For 3.5y, m Q (X) = 0 and F Q (X=3.5u) = 

—15 2 

7.1 x 10 watts/cm /y in accordance with the 
authors previously cited. Thus, equation 8 becomes 
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equation 9 . 

F nn „ (3.5u)~7.1 x 10 _15 watts lo"°' 4 m OBS (3 * 5p) (9) 

UBb ^ 

crn 


From the temperature and wavelength, the 
intensity B(A,T) of a black body radiator in units 
of power per unit area per steradian per unit wave- 
length interval may be calculated from the Planck 
function (equation 10) , 

C 1 

B ( A , T ) = - ^ (10) 

5 '"'2 

7T A D (e ~ -1) 

where and C 2 are constants,! is the wavelength 

-12 2 

and T is the temperature. C^=3.74l2 x 10 watts/cm 
and C 2 = 1.43879 cm deg (Allen, 1963). 

For the 3.5u wavelength, equation 10 becomes 
equation 11 . 


B (3. 5]i ,T) 


22.67 watts 

/ 4111 .. 2 

(e — -1) cm u ster. 


(ID 


From F_ c (3.5y) and B(3.5y,T), the solid angle 

UJdo 

subtended by the source in steradians (fi) may be 
computed by equation 12. 


Q = f qbs ( 3.5y)/B(3.5y ,T) 


(12) 


The angular diameter of the source (9) in milli- 
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seconds of arc may then be deduced via equation 13. 


-17 V2 

9 ms = ^/ 1 * 846x10 J 


(13) 


From the temperature and solid angle of the 
source, the apparent bolometric and visual magnitudes 
are next on the list of parameters to be found. 

The apparent bolometric magnitude (n' bol ) is found 
through equation 14, 

"W = - 2 - 5 lo 9 (F bol /F o bol> < 14 


where F q is the flux for zero apparent bolometric 

magnitude, F q = 2.52 x 10 ^watts/cm^ (Allen, 

1963) and F^^ is the computed total flux found 
from equation 15, 


F 

bol 



(15) 


where a is the Stefan-Boltzmann constant, 

{<7=5.70x10 watts/cm /deg 4 (Allen, 1963)} , T c is 
the color temperature in degrees Kelvin from 
equation 7 and Q is the solid angle from equation 12. 

The calculation of bolometric magnitude or lumin- 
osity from the derived color temperature assumes that 
the star radiates as a uniform disk, and limb darkening 
is ignored. The absolute value of the bolometric 
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magnitude may then be in error. However, the desire 
in the computations was to ascertain the shape of 
the luminosity curve as a function of phase and its 
relative range. 

The apparent visual magnitude, m , can also be 
determined from equations 14 and 10. The visual 
magnitudes were computed for an effective wavelength 
of 0.55y and the zero magnitude flux was taken to be 
F q (A=0.55u) = 3.92x10 12 watts/cm 2 /y (Johnson, 1965). 

Appendix IV lists the stars for which 1.04y and 
3 . 5y color indices exist as a function of visual 
phase. The first column tabulates the phase, the 
second column gives the color index, the third and 
fourth columns list the color temperature and angular 
diameter, columns 5 and 6 list the computed apparent 
bolometric and visual magnitudes, whereas column 7 
tabulates the mean observed visual magnitude. The 
data for column 7 were taken from the mean visual 
curves of Campbell (1955) . 

Figure 12 graphs the color temperature (T c ) , 
angular diameter (0) in milliseconds of arc, and 
apparent bolometric magnitude (m^Q^) as functions of 
phase. The data were obtained from appendix IV. 

R And, W And, R Cnc, S CrB, R Ser are omitted from 
the plot because the observational data were not 
complete over more than about one-half a cycle. Of 
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Figure 12c. Similar to figure 12a but for R Boo. 
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Figure 12d. Similar to figure 12a but for R Cas . 
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Figure 12e. Similar to figure 12a but for T Cas. 
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Figure 12f. Similar to figure 12a but for T Cep. 
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Figure 12g. Similar to figure 12a but for o Cet. 
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Figure 12i. Similar to figure 12a but for x Cyg. 
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Figure 12 j. Similar to figure 12a but for U Her. 
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Figure 12k. Similar to figure 12a but for R Leo. 
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Figure 12n. Similar to figure 12a but for U Ori. 
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Figure 12o. Similar to figure 12a but for R Tri. 
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Figure 12p. Similar to figure 12a but for R Vir. 
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course, the carbon rich stars are not plotted because 
no 1.04u data are available. 

There exists one point which should again be 
stressed at this time. The 1.04u and 3.5y magnitudes, 
from which the color indices were formed, represent 
an average behavior of the stellar variability. 

The color indices were formed from smoothed curves. 
During phases where either one or both of the 
magnitudes are rapidly changing, the color index is 
more likely to be in error. This error may be propa- 
gated through the computations and result in departures 
from a smooth curve for the temperatures or radii. 

It is extremely difficult to estimate the absolute 
accuracy of the results quoted in appendix IV, 
because the calculations are based on previously 
averaged data. As can be seen from prior figures 
(e.g. 4), the instantaneous brightness of a certain 
star does deviate from the adopted mean curves by up 
to + 20%. The mean adopted magnitudes hopefully 
represent the true average behavior of the star. 

The absolute accuracy of the color temperatures 
and angular diameters depends on a multitude of assumed 
values for various parameters. For example, if the 
assumption of a 10,000 °K color temperature for a Lyra 
is seriously in error, the color temperatures will be 
wrong. Other errors may be implicit: the absolute 
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calibration flux for zero magnitude at 3.5y, the 3 . 5y 
effective wavelength itself, the uniform disk 
approximation, as well as others. 

On a positive note, however, the relative un- 
certainties in the calculated parameters should be 
smaller than the absolute errors. One of the chief 
sources of uncertainty lies in the assumption that 
the 1.04y and 3 . 5y magnitudes represent the same 
optical depth in the stellar photosphere and are 
unchanged throughout the cycle. A second source of 
misgiving results from the averaging process used to 
determine the mean 1.04y and 3 . 5y magnitudes as 
functions of phase. 

Due to the above reasons, no error bars are plotted 
on the curves of figure 12. A conservative guess on 
the relative uncertainties of individual points is 
estimated to be about + 20%. 

Some of the individual stars of appendix IV and 
figure 12 merit further comment (other than the obvious 
ones) : 

T Cep: Although the shapes of the temperature and 
angular diameter curves are in fair agreement with 
other stars, the luminosity remains constant in direct 
contrast to the other variables. 

o Cet: Mira itself unfortunately does not possess the 
best quality 3.5y data. Its results should be viewed 
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with a more skeptical eye than most of the others. 

R Cyg and x Cyg : S type stars . Temperatures 

near minimum may be influenced by unknown bands 
affecting the 1 . 0 4y magnitudes (Lockwood, 1972). 

These temperatures must be referred to as equivalent 
temperatures . 

A few general comments are suggested by 
figure 12. The color temperature reaches a maximum 
near zero phase (maximum visual brightness) and attains 
a minimum value near phase 0.5, or approximately 
minimum visual magnitude. The angular diameters 
are largest near phase 0.3 and smallest near phase 
-0.1. The total luminosity maximum occurs near phase 
0.1 and minimum near phase -0.4. All of the above 
phases of extrema are in excellent agreement with 
Pettit and Nicholson (1933) , which appears to be 
the only other source of information for these 
quantities. A capsule description of Mira variables 
is that they are hottest when smallest and coolest 
when largest. 

Table XIV summarizes the average spectral type 
as well as the calculated quantities of temperature, 
angular diameter and luminosity extrema. The average 
spectral types are taken from table IV. Also included 
in table XIV are the average temperature and angular 
diameter, their absolute and percentage change, and 



Table XIV 

Summary of Calculated Parameters of the Miras 


Average Color Temperatures 


Star 

Spectral 

Type 

max 

min 

T 

C 

max 
0 K 

T c . 
min 

°K 

AT 

°K 

T 

°K 

AT 

T 

R Tri 

M4e + 

M8.2 

2850 

2390 

460 

2620 

.18 

R Boo 

M4 . 5e 

M8.0 

2750 

2480 

270 

2615 

.10 

R Vir 

H4.5e 

- 

2920 

- 

- 

- 

- 

o Cet 

M5.5e 

M9.0 

2790 

2320 

470 

2555 

.IS 

R Aql 

M6.5e 

M8.5 

2770 

2280 

490 

2525 

.19 

R Aur 

M6.5e 

M9.0 

2480 

1960 

520 

2220 

.23 

R One 

M6.5e 

- 

2650 

- 

- 

- 

- 

T Cep 

M6 . 5e 

M8.5 

2660 

2260 

400 

2460 

.16 

U Her 

M6.5e 

M9.0 

2510 

2010 

500 

2260 

.22 

x Oph 

MS . 5e 

M8. 5 

2680 

2440 

240 

2560 

.09 

0 Ori 

M6 . 5e 

M9.0 

2510 

2170 

340 

2340 

.15 

R Ser 

M6 . 5e 

- 

- 

2220 

- 

- 

- 

R Cas 

M7g 

M10 

2530 

1860 

670 

2195 

.31 

R Leo 

M7e 

M9 . 0 

2640 

2100 

540 

2370 

.23 

R LMi 

M7e 

M8. 8 

2590 

- 

- 

- 

- 


Luminosity 

Angular Diameter 


0 

max 

ms 

6 min 

ms 

A 0 
ms 

6 

ms 

A0 

9* 

m 

max 

Bol 

min 

Am 3ol 

7.2 

6.2 

1.0 

6,7 

.15 

3?5 

4?1 

0?6 

4.5 

3.7 

0.8 

4.1 

.20 

4.6 

5.2 

0.6 

3.9 

- 

- 

- 

- 

4 .6 

- 

~ 

40.2 

36.0 

4.2 

38.1 

.11 

-0.1 

0.7 

0.3 

17.1 

14.0 

3.1 

15.6 

.20 

2.0 

2.6 

0.6 

20.0 

14.5 

5.5 

17.3 

.32 

2.2 

2.9 

0.7 

- 

- 

- 

- 

- 

2.2 

- 

- 

27.5 

21.9 

5.6 

24.7 

.23 

1.3 

1.4 

0.1 

15.7 

13.8 

1.9 

14.8 

.13 

(2.4) 

3.3 

0.9 

17.1 

14.7 

2.4 

15.9 

.15 

2.0 

2.4 

0.4 

19.8 

14.5 

5.3 

17.2 

.31 

1.9 

2.9 

1.0 

- 

- 

- 

- 

- 

- 

3.8 

- 

32.5 

25.6 

6.9 

29.0 

.24 

1.0 

2.0 

1.0 

47.6 

34.9 

12.7 

41.2 

.31 

0.2 

0.9 

0.7 M 

17.4 

13.9 

3,5 

15.7 

.22 

2.2 

(2.8) 

(0.6) S 



Table XIV (continued) 


Color Temperatures 



Spectral 


c 

max 

c rain 

AT 

T 

AT 


Type 







Star 

max 

rain 

°K 

°K 

°K 

°K 

T 

T Ca.s 

M7. 5e 

M8.5 

2590 

2160 

430 

2375 

• 18 

R Cyg 

S3 , 9e 

- 

2460 

2060 

400 

2260 

.18 

X Cyg 

S4 , 7e 

M9 • 0 

2340 

1740 

600 

2040 

• 29 


Angular Diameter 

Luminosity 


e 

m§ a * 

0 . 
ms 

Ae 

IT1S 

9 

ms 

AG 

W 

ln Bol . 
max mm 

Al?l Bol 

20.7 

16.1 

4.6 

18.4 

.25 

1.9 

2.2 

0.3 

8.5 

4.7 

3.8 

6.6 

.57 

3-9 

5-3 

1.4 

44,4 

29.6 

00 

37.0 

• 40 

0-8 

1.8 

1.0 


u 

o 
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the range in total luminosity. The temperature (AT c ) , 
diameter (A0) and luminosity ranges are simply 

the differences between maximum and minimum, whereas 
the average temperature (T ) and average angular 
diameter (0) represent the arithmetic average of 
the extrema. 

Table XV lists the phases at which the color 
temperature, angular diameter, and luminosity reach 
their extreme values, as determined from figure 12. 
Also included in the last line are the mean values, 
the standard deviation of the mean and the number of 
stars used to compute the mean. 

Because some stars were radically different in 
the phases of the extrema of some of the parameters, 
they were omitted from the averages. U Her was not 
included in the angular diameter averages, o Cet 
was omitted at minimum 0, and, of course, T Cep was 
deleted from the bolometric magnitude phases because 
it appears to be constant over the entire cycle. 

From the results listed in table XIV, it comes 
possible to extend the color temperature-spectral 
type relation of table XII to lower temperatures and 
later spectral types. It must be emphasized, however, 
that the temperature versus spectral type function 
(as determined by non-variable or mildly variable M 
giants) may be inherently different from that found by 
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Table XV 

Phases at which Color Temperature, 
Angular Diameter, and Bolometric Magnitude 

Attain Extrema 


Star 

j. 

c 

max 

min 

max 

min 

Bol 

max 

'min 

R Aql 

0.00 

-0.48 

0.25 

-0.10 

0.07 

0.48 

R Aur 

-0.10 

0.40 

0.35 

-0.15 

0.12 

0.50 

R Boo 

0.10 

-0.42 

0.15 

-0.20 

0.12. 

-0.40 

R Cnc 

-0.05 

-- 

>0.25 

<0.15 

0.05 

— 

R Cas 

-0.05 

-0.42 

0.40 

-0.07 

0.10 

0.42 

T Cas 

0.00 

0.40 

0.27 

-0.07 

0.10 

0.45 

T Cep 

0.00 

0.50 

0.47 

-0.02 

Constant 

o Cet 

0.05 

-0.33 

0.25 

0.47 

(-0.02) 

0.15 

-0.40 

R Cyg 

-0.07 

0.50 

0.25 

-0.20 

0.12 

-0.32 

X Cyg 

-0.05 

0.50 

0.27 

-0.12 

0.10 

-0.45 

U Her 

0.00 

-0.35 

-0.30 

0.47 

(-0.03) 

>0.00 

-0.42 

R Leo 

0.00 

0.45 

0.35 

-0.05 

0.05 

-0.42 

R LMi 

0.00 

— 

0.25 

-0.10 

0.10 

— 

X Oph 

0.00 

0.38 

0.22 


0.07 

-0.47 

U Ori 

0.00 

-0.40 

0.17 

-0.22 

0.15 

-0.37 

R Ser 

— 

-0.45 

<0.2 5 

— 

— 

-0.45 

R Tri 

0.02 

-0.45 

0.50 

-0.05 

0.15 

-0.40 

R Vir 

0.15 

— 

0.30 

— 

-.20 

— 

Mean 

0.00 

-0.48 

0.30 

-0.11 

0-11 

-0.43 


+ 0.01 

+ 0.02 

+ .03 

+ .02 

+ .01 

+ .02 


N=17 

N=15 

N=15 

N=12 

N=15 

N=14 
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studying the Mira type variables. Table XVI lists 
the spectral type, color temperature and color index 
as deduced from the Miras in table XIV. The color 
temperatures for the Miras are a weighted average 
of the stars of each spectral type. Temperatures 
for spectral types M4 - M7 . 5 were determined from the 
average spectral type at maximum, while the tempera- 
tures for M8.0 - M10 were deduced from the average 
spectral type at minimum. 

-1 

The [1. 04u] - [3 . 5uJ color index versus T of 

c 

figure 13 is an extension and duplication of figure 11. 
The solid line is the assumed color index as a function 
of temperature for a 10,000 °K reference color tempera- 
ture which was defined to have zero index. The 
spectral types labeled to the left of the line 
represent the color index as determined from the normal 
M giant stars of table XII; the spectral types plotted 
to the right of the line are inferred from table XIV 
for the Miras. 

As is evident from figure 13, the temperature or 
color index agreement is best near spectral type M8 
and poorest for earlier types. As previously stated, 
interstellar extinction has been ignored. This 
reddening may affect the calculated temperatures but 
probably not enough to account for the observed 
discrepancy . 
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Color 

Table XVI 

Temperature and Color 

Index 

Versus Spectral Type for the Mira 

Variables 

Spectral 

Type 

Color Temperature 

Color Index 

M4-M4. 5 


2840 °K 

2 I ?0 8 

M5.5 


2790 

2.13 

M6. 5 


2610 

2.36 

M7 


2590 

2.40 

M7.5 


2590 

2.40 

M8.0 


2480 

2.56 

M8.2 


2390 

2.70 

M8.5 


2290 

2.90 

M8 . 8 


- 

- 

M9 . 0 


2110 

3.24 

M10 


1860 

3.87 



COLOR INDEX [l.04/i] - [ 3.5/i] 


135 



Figure 13. A comparison of the color index-spectral 
type relation versus color temperature as 
determined for the non-variable or mildly 
variable M giants and for the Miras. 



136 


A comparison between the effective temperatures 
for M giants, as listed in table XII, and the color 
temperatures for Miras (table XVI) , shows that the 
temperatures of Lee (19 70) , Johnson (1966) and Smak 
(1964) bracket those of the present study. Jones 
(1968) produced a color temperature-spectral type 
relation for M giant stars from narrow band photometry 

O 

in the 4700-11200 A region. When his color tempera- 
tures are compared to those deduced from the non- 
variable and variable stars of this study, the 
agreement is fair. The two extremely different 
wavelength regions do not yield discordant color 
temperatures . 
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VIII 

VISUAL MAGNITUDES AND BLANKETING 

In addition to the previous calculations of 
color temperature, angular diameter, and bolometric 
magnitude, a comparison between the calculated and 
the observed visual magnitudes allows a deter- 
mination of visual absorption which is caused by 
lines or molecular bands. 

The first step in this procedure is to deduce 
the visual blanketing for the non-variable or mildly 
variable M giant standard stars of table XI. From 
the adopted color index-color temperature relation, 
temperatures dependent on spectral type may be 
calculated in a manner similar to that previously 
described. The earliest spectral types possess color 
indices less than 2 1 ?0 . For these color indices , the 
linear approximation between color index and T c 
is no longer valid and it was not used. The color 
temperatures for these "early" sources were calculated 
by interpolation of table XIII. Once the color 
temperature is known, the solid angle subtended by 
the source, the uniform disk angular diameter, and 
the visual magnitude are computable. Of course, the 
implicit assumption in these calculations is that 'the 
stellar photosphere is radiating as a black body wftich 
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may be characterized by one temperature. This unique 
temperature is assumed to be applicable to all three 
(0.55, 1.04 and 3 . 5y ) wavelengths. 

Table XVII again lists the standard M giants, 

their spectral type, color index (C.I.), color 

temperature (T ) , angular diameter (0), calculated 

and observed visual magnitudes (m . , m , ) and 

v calc v obs 

the difference (A) between them. 

As is evident from table XVII, MO stars appear 
brighter than their calculated visual magnitude by 
about 0^3. This difference decreases, changes 
sign about M3 and then grows larger as the spectral 
type advances . 

This apparent excess of MO-M2 stars can be 
attributed to a variety of causes . One source of 
discrepancy may be in the adopted 3.5y effective 
wavelength and its absolute calibration. However, 
for the purposes of the discussion which follows, 
the disagreement between the observed and calculated . 
visual magnitudes near spectral types MO and Ml is 
considered remarkably small and will not alter the 
ensuing results. The main point to be stressed is 
that as the temperature decreases, the amount of 
blanketing at visual wavelengths increases. 

The next step in the visual calculations is to 
determine the blanketing in the same manner for the 



Table XVII 

Calculated Parameters for M Giant Standard Stars 



Spectral 

Color 

T 

C 

0 

m 

V 

m 

V 


Star 

Type 

Index 

°K 

ms 

Calc 

Obs 

A 

3 And 

MO 

III 

-.m, . 
1.64 

3289 

16.6 

2?41 

2™05 

-o™36 

M UMa 

MO 

III 

1.60 

3333 

10.2 

3.37 

3.04 

-0.33 

75 Leo 

MO 

Illb 

1.60 

3333 

3.7 

5.57 

5.18 

-0.39 

a Vul 

MO 

III 

1.58 

3356 

5.3 

4.74 

4.45 

-0.29 

Y Eri 

Ml 

III 

1.51 

3390 

9.9 

3.28 

2.95 

-0.33 

X Peg 

M2 

III 

1.57 

3367 

5.0 

4.83 

4.80 

-0.03 

a Cet 

M2 

III 

1.63 

3300 

15.0 

2.61 

2.53 

-0.08 

ir Leo 

M2 

III 

1.78 

3125 

6.2 

5.02 

4.71 

-0.31 

83 UMa 

M2 

Illab 

1.66 

3257 

5.8 

4.77 

4.73 

-0.04 

3 Peg 

M2 

Ilb-IIIa 

1.80 

3106 

21.0 

2.41 

2.53 

0.12 

55 Peg 

M2 

III 

1.59 

3344 

5.2 

4 .79 

4 . 50 

-0.29 

p Gem 

M3 

Ilia 

1.76 

3155 

17.3 

2.69 

2.85 

0.16 

p Per 

M4 

Ilia 

1.84 

3067 

18.6 

2.79 

3.39 

0.60 



Table XVII (continued) 


Star 

Spectral 

Type 

Color 

Index 

T 

c 

°K 

e 

ms 

*n 

V 

Calc 

V 

Obs 

A 

w Vir 

M4-4.5 III 

1.79 

3125 

7.9 

4.48 

5.28R 

0.80 

HR 5299 

M4-4.5 III 

1.87 

3040 

8.9 

4.47 

5.28 

0.81 

VY Leo 

M5 III 

1.96 

2941 

11.3 

4.28 

5.78 

1.50 

RR UMi 

M5 III 

1.76 

3155 

10.8 

3.71 

4.59 

0.88 

4 

x Ser 

M5 Ilb-IIIa 

2.23 

2717 

13.8 

4.63 

— 


OP Her 

M5 Ilb-IIIa 

1.93 

2976 

7.3 

5.10 

6.32 

1.22 

R Lyr 

M5 III 

1.96 

2941 

20.6 

2.98 

4. OH 

1.02 

W Cyg 

M4e~M6 

2.17 

2762 

16.7 

4.06 

5.1H 

1.04 

RZ Ari 

M6 III 

2.02 

2890 

12.9 

4.16 

5.93 

1.77 

RS Cnc 

M6 Ilia 

2.16 

2770 

18.8 

3.77 

5.95 

2.18 

X Her 

M6e 

2.07 

2841 

15.9 

3.88 

— 


g Her 

M6 III 

2.07 

2 841 

21.2 

3.25 

5.05 

1.80 

CH Cyg 

M6 

2.48 

2532 

15.0 

5.22 



SW Vir 

M7 III 

2.23 

2717 

19.5 

3.89 

- 

— 

RX Boo 

M8 : 

2.64 

2433 

24.9 

4.58 

- 

- 
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Mira variables. Appendix IV lists the calculated 
and observed visual magnitudes as a function of phase 
(or alternatively, of temperature) for the Miras. 

As stated previously, the observed brightness is 
the mean visual magnitude as determined by Campbell 
(1955) from the AAVSO observations. The effective 
wavelength for this data (based on the human eye) 
may be somewhat different than the 0.55u used in the 
calculations. This difference may introduce uncer- 
tainties into the absolute amount of visual absorption. 
However, with the lack of simultaneous UBV photometry 
or mean photometric light curves for the entire cycle, 
there is no alternative. 

In order to compare similar points in the stellar 
cycle, the times or phases of maximum and minimum 
temperature were chosen as reference marks. These 
times were then assumed to represent earliest and 
latest spectral types. The differences between the 
observed and calculated visual magnitudes were 
tabulated for these reference points. 

As an alternative method of computation, the 
reference marks were selected to be the phases of 
visual observed extrema (e.g., zero phase for 
maximum visual magnitude and near phase +0.50 for 
minimum) . The difference in visual magnitude was also 
tabulated for these points. 
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Table XVIII summarizes the results of these two 
different methods. Column 1 lists the spectral type, 
column 2 names the star whose average spectral type 
at maximum or minimum fulfills column 1, column 3 
lists the phase of the temperature extrema (maximum 
T c for spectral types M4 - M7.5; minimum T c for later 
spectral types) , column 4 tabulates the difference 
between observed and calculated visual magnitudes. 

Column 5 lists the phase of visual observed extrema 
and column 6 tabulates the magnitude difference for 
this alternative method. Due to the visual influence 
of its companion, X Oph was omitted at minimum from 
the calculations. 

As evidenced by table XVIII , there is really 
very little difference in these two alternative 
methods of computation. For this reason, both results 
were used to determine a mean visual absorption as a 
function of spectral type for the Miras. 

Smak (1964) attempted to ascertain the amount of 
visual blanketing by TiO bands as a function of spectral 
type or effective temperature for M giants. His 
results for blanketing and his effective temperatures 
are listed in table XIX. Included in table XIX, and 
as determined by this study, are the blanketing and 
color temperature (table XVII) for the non-variable 
standard stars. These same quantities for the Miras, 



Table XVIII 
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Visual 

Blanketing Summary 

for Miras 


Spectral 

Type 

Star 

Phase of 

T extrema 
c 

A 

Phase of 
Visual 
Extrema 

A 

M4 

R Tri 

0.02 

0™7 

0.00 

0^5 

M4.5 

R BOO 

0.10 

0.8 

0.00 

0.1 


R Vir 

0.15 

0.9 

0.00 

0.0 

M5. 5 

o Cet 

0.05 

1.6 

0.00 

1.4 

M6.5 

R Aql 

0.00 

1.9 

0.00 

1.9 


R Aur 

-0.10 

3.8 

0.00 

2.2 


R Cnc 

-0.05 

2.2 

0.00 

1.9 


T Cep 

0.00 

2.0 

0.00 

2.0 


U Her 

0.00 

2.2 

0.00 

2.2 


X Oph 

0.00 

2.3 

0.00 

2 . 3 


U Ori 

0.00 

1.5 

0.00 

1.5 

M7 

R Cas 

-0.05 

2.7 

0.00 

2 . 4 


R Leo 

0 . 00 

3.0 

0.00 

3.0 


R LMi 

0.00 

2.1 

0.00 

2.1 

M7. 5 

T Cas 

0.00 

3.0 

0.00 

3.0 

M8 

R Boo 

-0.42 

4.0 

0,50 

4.2 

M8 . 2 

R Tri 

-0.45 

4.0 

-0.47 

4.0 

M8 • 5 

R Aql 

-0.42 

5.3 

-0.45 

5.2 


T Cep 

0.50 

5.1 

0.47 

5.1 


T Cas 

0.40 

5.4 

0.45 

5.6 
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Table XVIII (continued) 


Phase of 


Spectral 

Type 

Star 

Phase of 

T extrema 
c 

A 

Visual 

Extrema 

A 

M9 .0 

R Aur 

0.40 

6.0 

0.50 

6.7 


o Cet 

-0.32 

4.9 

r- 

« 

o 

1 

4.9 


U Her 

-0.35 

3.9 

l 

o 

• 

O 

4.3 


R Leo 

+0.4 5 

4.3 

-0.42 

4.7 


U Ori 

-0.40 

4.8 

-0.37 

4.8 


R Cas -0.42 


M10 


4.4 


-0.40 


4.5 
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Table XIX 

Visual Blanketing as a Function of Spectral Type 


Spectra 

Type 

Smak (1964) 

T <5 T , 0 

e i 

Non- 

T 

c 

THIS 

-Variables 

A 

STUDY 

Miras 

T A 

c 

MO 

3340 

0™00 

3320 

-0™34 



Ml 

3200 

0.22 

3440 

-0.33 



M2 

3090 

0.45 

3240 

-0.10 



M3 

2980 

0.77 

3150 

0.16 



M4 

2850 

1.13 

3070 

0.74 

2850 

0?6 

M4 . 5 





2840 

0.4 

M5 

2710 

1.59 

2940 

1. 15 



M5 . 5 





2790 

1.5 

M6 

2500 

2.27 

2840 

1.92 



M6.5 





2610 

2.1 

M7 

2250 

3.23 

2710 

- 

2590 

2 . 6 

M7 . 5 





2590 

3.0 

M8 

2000 

4.33 

2430 

- 

2480 

4.1 

M8 . 2 





2390 

4.0 

M8.5 





2290 

5.3 

M9 





2110 

4.9 

M10 





1860 

4.5 
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as summarized in tables XVI and XVIII, are also listed. 

As may be deduced from table XIX, this color 
index-color temperature method agrees quite well 
with Smak's (1964) values. The difference in STiO 
near spectral type MO cannot be simply resolved by an 
additive constant but is deeply embedded in the 
absolute fluxes, color temperatures and other 
assumptions previously quoted and cannot be further 
resolved at this time . However , to obtain agreement 
between the two summaries , the addition of about 
0.3 to 0.4 magnitudes to column 5 or, as an alter- 
native, the subtraction of that same amount from 
column 3, would bring the results into excellent 
agreement. It deserves to be mentioned that Smak 
defined his STiO correction to be zero at spectral 
type MO . 

In addition to the 6TiO absorption at temperature 
or visual magnitude extrema, the visual blanketing as 
a function of phase, or probably more significantly 
as a function of color temperature, may be computed 
throughout the entire cycle of variation. Figure 14 
depicts the <STiO absorption (A=m^ ca i c -in v Q bs^ versus 
color temperature for the M type Miras which possess 
complete cycle parameters. The data for figure 14 
was taken from appendix IV. The binary, X Oph, is 


omitted . 



VISUAL ABSORPTION A ( MAG ) 



2000 2200 2400 2600 2800 3000 


COLOR TEMPERATURE T c 

Figure 14a -Visual absorption in magnitudes versus 
color temperature throughout the 
full cycle of variability for the 
earliest M stars of this study. 



2000 2200 2400 2600 2800 3000 
COLOR TEMPERATURE T c 

Figure 14b. Similar to figure 14a but for M stars 
of later spectral type. 



VISUAL ABSORPTION A (MAG) 
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COLOR TEMPERATURE (T c ) 

Figure 14c. Similar to figure 14b. 



VISUAL ABSORPTION A (MAG) 
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2000 2200 2400 2600 2800 3000 


COLOR TEMPERATURE <T C ) 


Figure 14d. Similar to figure 14b. 


VISUAL ABSORPTION A (MAG) 
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Two S stars ,x Cyg and R Cyg, have been included 
in figure 14 for comparative purposes. As previously 
mentioned, near minimum light, the 1.04P magnitudes 
are influenced by unidentified absorption bands. 

The color temperatures, as deduced from the color 
indices, may then be cooler than appropriate for 
the photosphere. Other quantities calculated from 
this temperature would also be in error. 

Lockwood and Wing (1971) give no equivalent M 
spectral types for R Cyg but cite equivalent types 
M2-M9.0 for the mild S star, x Cyg. For x Cyg, the 
very early spectral type is consistent with the rather 
small amount (-.374) of visual absorption at zero 
phase. The 4™7 of absorption at minimum temperature 
is also in agreement with the M9 . 0 spectral type. 

The visual blanketing-color temperature relation 
throughout the complete cycle appears similar to that 
of the M stars. 

As is evidenced by figure 14, the visual blanketing 
is not a unique function of color temperature. There 
appear to be two distinct values of absorption for a 
given temperature, depending on whether the star is on 
the ascending or descending branch of the light curve. 
The amount of blanketing executes loops in the blanket- 
ing-temperature diagram akin to a hysteresis curve. 

This behavior has been previously noticed by Spinrad 
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and Wing (1969) in their plots of molecular band 
strength versus temperature. 

In figure 14 , not all of these loops circle 
in the same sense. The clockwise loops imply 
increasing absorption as the color temperature 
decreases on the descending portion of the light 
curve. As the color temperature reaches a minimum 
value and begins to rise, the blanketing does not 
recover as quickly, possibly implying a further 
formation of, or increase in, the amount of TiO. 

This amount of TiO is then not a unique function 
of temperature, but depends on the stellar history. 

To the contrary of these clockwise loops , a 
few stars cycle in the opposite sense. It is 
perhaps significant that these exceptions tend to 
be of earlier spectral type. Spinrad and Wing (1969) 
noticed these opposite loops also. 

The exact processes contributing to these loops, 
as well as the significance of the clockwise or 
counter-clockwise behavior, are beyond the scope 
of this present work. 
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IX 

RADIAL VELOCITIES 

Another quantity which may be computed in order 
to compare with observational data is the radial 
velocity of the stellar photosphere. Scott (1942) 
assumed that the bright hydrogen emission lines of 
Mira variables originated at about the same level 
as the photosphere. His calculations, based on this 
assumption, provided an acceptable explanation for 
the luminosity variation of Mira. 

Unfortunately, the emission lines are present 
only for about one-half cycle centered on maximum 
light, and, therefore, cannot be used to trace the 
photospheric radial velocity throughout the entire 
cycle. Merrill and Burwell (1930) published radial 
velocity curves of emission lines for selected 
Miras. The stars which are in common with this 
project are: T Cas , R Ser , U Ori , R Leo, x Cyg , 

R Cyg , and R Cas. T Cas, in contrast to most other 
stars , had either a constant velocity or a smaller 
range. The trait of the velocity curve in common 
among the Miras appears to be a relatively broad 
region of maximum velocity of approach, sometime 
after zero phase. This expansion velocity seems to 
maximize somewhere near phase 0.15. Joy (1926) 
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produced a mean curve for the emission line velocity 
for Mira itself in which the expansion velocity 
maximizes near phase one-quarter, minimizes near 0.5 
phase and has secondary minima and maxima near 
phase 0.9 and 0.75, respectively. 

If the assumption of emission lines originating 
near the photosphere is valid, the radial velocities, 
as calculated from the angular diameters in appendix 
IV, should approximate the emission line velocities 
where comparisons can be made. Due to their apparent 
brightness and data coverage, three representative 
stars (R Leo, R Aql, and T Cas) were selected for 
comparison purposes . The radial velocities were 
computed from equation 16 to represent the shape of 
the velocity curve, not its absolute scale. 

e(*) = 0 (*+o.05)-e (4>) + [e (40 -e (*-o.os) ] 

= e (4>+o.05)-e (*-o.05) (16) 

% 

where 0 represents the radial velocity in arbitrary 
units and * is the phase. This equation was utilized 
in an attempt to average or smooth out the fluctuations 
and errors in computation of temperature and angular 
diameters. Figure 15 graphs this radial velocity as 
a function of phase for R Leo, R Aql, and T Cas. 
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Figure 15. Computed radial velocity versus phase 
for three representative stars R Leo, 
R Aql , and T Cas . 
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A positive 0 implies expansion. 

The general shape of the radial velocity curves 
is emphasized, not the detailed structure. The 
angular size absolute error estimates are already 
quite large, although the relative errors should be 
smaller. From figure 15, it is apparent that the 
outward velocity reaches maximum near phase +0.1, 
while the inward velocity attains a minimum near 
phase +0.45. Secondary velocity extrema occur near 
phase -0.4 and -0.3, respectively. The maximum 
diameter takes place near phase 0.3 where the outward 
velocity reaches zero, while the minimum diameter 
occurs near phase -0.1. Both these results are in 
agreement with table XV -- as they should be. They 
are also in reasonable agreement with Merrill and 
Burwell (1930). 

It is unfortunate that the 3 . 5y data for o Cet 
are not more complete or trustworthy so that a 
better radial velocity comparison may be made on the 
same star. The general shape of the radial velocity 
curve of R Leo does appear similar to the hydrogen 
emission line velocity curve of Mira as determined by 
Joy (1926). A shift of R Leo's velocities by about 
0.1 cycle toward larger phase would bring the two 
curves into excellent agreement. The two curves show 
the primary maximum of outward velocity, the primary 
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minimum, and suggestions for the secondary extrema. 

The reason for the apparent displacement of the 
computed radial velocities of figure 15, with respect 
to the observational velocities, is not clear at this 
time. The calculated results are not in good agreement 
with Joy's (1926) emission lines, nor Scott's (1942) 
corrected values. The deduced velocities of figure 15 
do not mirror the behavior of the absorption lines. 

On the plus side, the velocity curves of figure 
are in good agreement with that deduced for Mira by 
Pettit and Nicholson (1933). From their graphs, 
the maximum expansion velocity happens near phase 
+0.1 and has a broad minimum centered near phase 


-0.4. 
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X 

CARBON-RICH MIRAS 

Although the Mira type and semi-regular 
variables of the C spectral designation have not 
been included in any of the previous computations, 
the characteristic behavior deserves to be mentioned 
at this time. 

The carbon stars observed as part of this project 
are listed in table II, but are enumerated again as: 

T Cnc, WZ Cas , S Cep, V CrB , U Cyg , RS Cyg , RY Dra , 

R Lep, W Ori , and SS Vir. Their period, spectral 
type, visual magnitude and 3.5y range were also 
summarized in table II. 

Unfortunately, the 3 . 5pi data are either 
incomplete or indecisive for T Cnc, RY Dra, R Lep, 

W Ori and SS Vir. T Cnc, RY Dra and W Ori are semi- 
regular variables of relatively small visual 
amplitude. Their 3.5U observations could be 
adequately described in each case by a constant 
magnitude (+0™26, -O^Sl, -0^89, respectively) with 
about a 10% uncertainty. The 3 . 5y data are in- 
sufficient on R Lep to allow any statement. Visual 
data were available for SS Vir, although the 3.5y 
points are of limited phase span. 
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Another semi-regular variable, RS Cyg , had 
visual observations available. However, even though 
the 3 . 5y coverage is well distributed over the entire 
cycle, it demonstrated no alterations beyond the 
observational uncertainties and its 3.5u magnitude 
is adopted as +0.94 + O^l. To the contrary, another 
small amplitude semi-regular variable, WZ Cas , 
showed significant 3 . Sp. changes, but no visual data 
for phase computation were available from the AAVSO. 

On a more positive note, three carbon stars, 
namely S Cep, V CrB and U Cyg, possess good 3 . 5y 
data coverage as well as AAVSO visual curves. Their 
3.5|i measurements were plotted in figure 6 and the 
adopted mean curves were listed in table V. 

In contrast to the M stars, no 1 . 0 4 p measurements 
throughout the variability cycle exist for carbon 
stars; therefore, no color temperatures, angular 
sizes nor luminosities may be similarly computed. 

These parameters must wait until observations which 
represent a true photospheric continuum level become 
available . 

Mendoza and Johnson (1965) have published 
UBVRIJKLM observations of carbon stars and concluded 
that the observed energy distributions "differ 
relatively little from those of ordinary K and M 
giants for wavelengths longer than that of the V 



filter (0 . 55y ) but that they differ greatly for the 
shorter wavelengths". They imply the effective 
temperatures of carbon stars are similar to the 
M giants. 

Gordon (1968) developed a spectral classification 
method for carbon stars and Richer (1971) deduced 
their color temperatures on the basis of the (I-L) 
index of Johnson (1966). Fujita et al . , (1965) 

published color temperatures for carbon stars and 
graphed a conversion from C to M spectral type. 

Table XX lists their conversion from C to M type, 
their color temperature, the color temperature of 
Richer (1971) and the color temperatures for M 
stars deduced from this study. As may be seen from 
the table, the agreement among the various methods 
of color temperature determination seems fair to good. 
If the conversion from C to M classification of 
Fujita, et al. , (1965) is valid, table XX supports 

Mendoza and Johnson's (1965) contention that the 
temperatures of carbon stars are similar to M stars. 

Other quantities existing at the present are 
the visual and 3 . 5y data which then may be compared 
for M, S, and C stars to ascertain similarities or 
differences. From the 3 . 5y phase curves of figure 6, 
there seems to be an indication that carbon stars may 
reach their 3.5y maximum somewhat nearer to zero phase 



Table XX 



Intercomparison of 

Carbon Stars and M 

Stars 


Spectral 

Type 

Color 

Temperature 

°K 

M Equivalent 
Spectral 

Type 

* 

Color 

Temperature °K 
This Study 

non- 
variables Miras 

CO 


4400 




Cl 


(3800) 




C2 

~ 32 50 

(3200) 

~MO 

3320 

— 

C3 

3050 

2700 

M2 

3240 

— 

C4 

2900 

2600 

M3. 5 

~ 3100 

— 

C5 

2750 

2450 

M5 

2940 

-2820 

C6 

2600 

2300 

M6 

2840 

-2700 

C 7 

2500 

2200 

M7 

2710 

2590 

C8 

2350 

2000 

M8 

2430 

2480 

C9 

2200 

1800 

M9 

— 

2110 

* T'irp-ta , 

e± al . , 1965. 





* _32±nher 

(1971) . 
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than do the M and S stars. The C type stars appear 
to minimize about the same phase as M stars but 
earlier than the S stars. Eight M stars exist for 
which a mean 3.5y range may be computed. For this 
group, the mean visual range is 5™3 + 0^2 (standard 
deviation of the mean) whereas the observed mean 
3 . 5y range is 0^5 + 0™1. For the four S stars, the 
mean visual range is 7 1 ?! + 0™4 and at 3 . 5y it is 
0™9 + 0^2 . For the carbon stars, excluding the 
constant RS Cyg, the mean visual range is 3^2 + 0™3 
while the 3 . 5y range is O 1 !^ +_ 0 1 ?!. Thus, the mean 
range of 3.5y variability appears very similar, 
although C and S type stars may have a somewhat 
larger range than do the M stars. The smaller 
visual range of the carbon stars is probably 
reflected in the excess TiO absorption of the M and 
S stars and is not likely due to inherently different 
temperatures . 

From the available data of this study, in agree- 
ment with Mendoza and Johnson (1965) , the carbon 
stars do not differ significantly from the M and S 
type variables. However, there is a suggestion, 
supported by only a few stars, that the carbon stars 
possess a somewhat greater range in 3.5y magnitude and 
tend to reach 3.5y maximum nearer zero phase than do 


the M type stars. 
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XI 

COMPARISONS WITH OTHER OBSERVATIONAL RESULTS 

Theoretical studies of long period variable stars 
are not only sadly lacking in number, but the pre- 
dictions which do exist do not lend themselves to a 
comparison of observable parameters. Therefore, 
this section will compare the computations of this 
present study with other observational results. 
Comparison of color temperature, effective temperature, 
visual absorption and radial velocity have been 
previously mentioned. 

The number of cool giants whose angular diameters 
have been measured is sparse. For the non-variable 
standards of table XVII, only (3 And, a Cet and 8 Peg 
have angular diameters determined by various inter- 
ferometric methods (Pease, 1931; Brown, 1968; 

Gezari, Labeyrie and Stachnik, 1972). Table XXI 
summarizes these observations as well as the uniform 
disk diameters calculated in this project and taken 


from table XVII. 
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Table XXI 

Summary of Angular Diameter Measurements 



Angular 

Diameter (ms) 

Star 

This study 

previous work 

8 And 

16.6 

10.8* 

a Cet 

15.0 

9.4 - 11.5* 

6 Peg 

21.0 

21:, 16 / 


* Brown (1968) (as extracted from the Carnegie Yearbooks) . 
: Pease (1931) . 

7 * Gezari, et al. , 1972. 

As table XXI indicates, the agreement is satisfactory 
only for 6 Peg, although the more recent measurements 
by speckle interferometry (Gezari et al., 1972) have 
reduced its size. 

For the Mira variables, only o Cet itself has 
been measured with the 20-foot Michelson interfero- 
meter. Its angular diameter was 47 ms (Pease, 

1931) , as compared to 36 - 40 ms computed by this 
study. Mira has also been measured by the speckle 
interferometry methods of Bonneau and Labeyrie (1973). 

They reported, at visual phase 0.09, an angular 

O 

diameter of 51 + 5 ms at a wavelength of 7500 A with 
the diameter increasing toward bluer wavelengths. 
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At 1.04y, they set an upper limit of 50 ms. At 
visual phase 0.38, they also reported a 55 + 5 ms 

O 

diameter at 7500 A. The diameters inferred by this 
project at phases 0.09 and 0.38 are 38 and 37 ms t 
respectively. The angular size of R Leo has been 
measured at infrared wavelengths by lunar occultation 
(Becklin, 1972). Its reported extent of 50 ms 
compares well with the 35-47 ms range of this study. 

As previously stated, the absolute errors of the 
computed angular diameters of this study are probably 
quite large (~30%) , although the relative errors should 
be smaller. Even though the agreement is far from 
excellent, it seems that the angular diameters deduced 
from the [l.Q4y] - [3.5y] color index and the 3.5y 
absolute calibration are not inconsistent with the 
diameters measured by other means. There is an 
indication (Bonneau and Labeyrie, 1973) that the 
diameters decrease at longer wavelengths. This 
reduction would bring the results for the Miras 
into better agreement but would compound the 
differences for the standards of table XXI. 

In a project similar to the present research, 
one of the earliest and, yet, most often quoted studies 
of LPV stars is the work of Pettit and Nicholson (1933) . 
Working with a thermocouple at the Newtonian focus 
of the Mt. Wilson 100-inch telescope from 1921-1927, 
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they measured 11 long period variables. From their 
observations of radiometric magnitudes and water 
cell absorption, they determined temperatures, angular 
diameters and luminosities as functions of visual 
phase . 

Comparison can be made of 9 stars (R Aql , R Cnc, 
o Cet, R Cyg,x Cyg , R Leo, R LMi , X Oph, R Tri) which 
the two projects held in common. Table XXII 
summarizes the computed temperature, diameter and 
radiometric magnitude range of Pettit and Nicholson 
(1933) and the similar quantities from this work. 

As is evident from the table, there exists some 
discrepancy as to the temperature and diameter 
extrema, although, in some cases the agreement is 
excellent. Differences in absolute calibration, 
but, particularly differences in methods of tempera- 
ture determination could account for the disagreement. 

Despite some of the individual star dis- 
crepancies which exist between this present study 
and the earlier work of Pettit and Nicholson (1933) , 
several corroborating generalizations concerning the 
average behavior of Mira variables can be deduced. 
Table XXIII tabulates the mean characteristics as 
determined by the two projects. The errors listed 
represent the standard deviation of the mean. Some 
of Pettit and Nicholson's results are quoted without 



Table XXII 

The Measurements of Pettit and Nicholson 
Compared to the Present Study 
Pettit and Nicholson This Study 


Star 

Temperature 
Range °K 

Angular 

Diameter 

ms 

Radiometric 

Magnitude 

Color 

Temperature 

0 

ms 

^bol 

R Aql 

2320-1920 

19-28 

1. 8-2.7 

2770-2280 

14-17 

2. 0-2. 6 

R Cnc 

2320- (1700) 

16- (31) 

1.9-2. 9 

2650- 


2.2- 

o Cet 

2640-1920 

38-55 

-0.3+0. 8 

2790-2380 

36-40 

-0.1-0. 7 

R Cyg 

2790-(2010) 

5.4- (9.7) 

3. 6-5.0 

2460-2060 

4.7-8. 5 

3. 9-5. 3 

X Cyg 

2240-1640 

34-57 

0.4-1. 6 

2340-1740 

30-44 

0. 8-1. 8 

R Leo 

2170-1800 

50-63 

0 . 1 - 1.0 

2640-2100 

35-48 

0.2-0. 9 

R LMi 

2270- (1780) 

18-24 

2. 0-3.1 

2590- 

14-17 

2.2- 

X Oph 

2390-1890 

19-28 

1.6-2. 3 

2680-2440 

15-17 

2. 0-2. 4 

R Tri 

3270-2460 

4. 3-7.8 

3.4-4 .3 

2850-2390 

6. 2-7. 2 

3. 5-4.1 
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Table XXIII 

Mean Variable Characteristics 


Pettit and 


Quantity 

This Study 

Nicholson 

Visual Range 

c m . . _m_ 
5 . 4 + 0 . 3 

5?9 

Luminosity Range 

0^7+0 . 1 

0™9 

AT/T 

0.19+0.02 

0.26 

Ae/'e 

0.25+0.03 

0.37 

Phase of Temperature 

Maximum 

0.00+0.01 

-0.01+0.01 

Minimum 

-0.48+0.02 

0.49+0.02 

Phase of Angular Diameter 

Maximum 

0.30+0.03 

0.28+0.02 

Minimum 

— 0 . 11+0 . 02 

-0.03+0.01 

Phase of Luminosity 

Maximum 

0.11+0.01 

0.14 

Minimum 

-0.43+0.02 

-0.45 
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error estimates because no uncertainty was published. 
Due to the smaller number of objects, the errors 
are probably at least as large as those calculated 
from this study. The Pettit and Nicholson data for 
the temperatures and diameters of 6 stars were taken 
from their table III, graphed, and statistically 
treated in a manner similar to the results of the 
present work. 

As is evident by table XXIII, the agreement 
between the two projects is astonishingly excellent. 
Although slight differences exist, the present 
study, with its improved infrared equipment, confirms 
the results of 50 years ago. 

Again, as evidenced by table XXIII, the maximum 
luminosity occurs near phase +0.1 and lags the maxi- 
mum compression by about one-quarter cycle. The 
number of stars is not sufficient to ascertain this 
phase lag of luminosity with respect to minimum radius 
as a function of composition, spectral type or 
visual range. 

The physical explanation for the lag of luminosity 
behind visual brightness can be explained as follows : 
the visual energy output is dependent on blanketing 
but is a very strongly dependent function of tempera- 
ture (e.g., exponential), whereas the luminosity is 
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dependent on the fourth power of the temperature. 

As the phase increases from zero, the temperature 

declines but the photosphere expands; and, for a 

short period of time (~0.1 cycle) , the solid angle 

4 , 

is increasing more rapidly than T is declining. 
This phase lag is also apparent at the 3.5y wave- 
length where the observed energy output is approxi- 

2 

mately proportional to R T. 
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XII 

MODELS FOR RED GIANTS 

Unfortunately , there has been a scarcity of 
theoretical models which attempt to explain the 
physical properties or the basis for the rather 
extensive observational data. Admittedly, the 
observations which have been gathered on Long 
Period Variable Stars are of a descriptive or catalog 
nature. Some models have emerged, but the predicted 
quantities may not be directly observable, and 
confirmation is difficult or non-existent. 

The dynamical models for LPV stars require 
static models of stellar structure for non-variable 
red giants as a springboard from which pulsation theory 
may be applied. Only in recent years have reasonable 
model atmospheres been constructed which take into 
account all of the recent infrared observations of 
these cool giants. The photospheric temperatures 
typically range from 2000 to 3500 °K and imply a 
Planckian energy distribution peaking near ly. 

Only since the early 1960 's have methods become 
available which allow the detailed structure of the 
cool stellar atmosphere to be explored. 

Various authors have provided important con- 
tributions to the theory of model atmospheres for 
red giants. Review articles by Munch (1960) on 
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stellar atmospheres; Keenan (1960) on physical 
properties; Merrill (1960) on spectra; Weymann (1963) 
on mass loss; Smak (1966) on general traits; Johnson 
(1966) on infrared measurements, temperatures and 
bolometric corrections; Spinrad and Wing (1969) on 
infrared spectra; Vardya (1970) on late type 
atmospheres; and Neugebauer, Becklin and Hyland 
(1971) on sources of infrared radiation, together 
provide the most recent information and theory 
regarding red giant characteristics. The interior 
structure of red giants has been calculated by Harm 
and Schwarzschild (1966) , while Iben (1967) has 
studied their evolutionary status. This preceding 
list is by no means complete and exhaustive, but 
it is included only as a starting point (or perhaps 
a terminating point) for the study of red giants. 

The physical state of a one solar mass red 
giant’s interior may be described as an inert, 
hydrogen deficient, isothermal core at a temperature 

of about 50 million degrees, a pressure of about 

21 2 -3 
10 dynes/cm , a radius of about 10 of the 

stellar radius and a mass of about one-fourth of 

the stellar mass (Harm and Schwarzschild, 1966). 

The energy generation arises from hydrogen burning 

in a very narrow shell surrounding the inert core. 

The energy transport is by radiation until convection 
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sets in at about 1/50 of the stellar radius. At 

the photosphere of this one solar mass red giant, 

typical values are: radius -2 00 R^, temperature 

2 - 1 0 ^ 

-3000 °K, pressure -20 dynes/cm , density -10 gms/cm , 

2 

and surface gravity -0.3 cm/sec (Keeley, 1970a). 

Because this present work deals primarily with 
stellar variability, the red giant models have only 
been briefly discussed in order to provide background 
information as to their physical nature. 
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XIII 

MODELS FOR LONG PERIOD VARIABLE STARS 

It is now generally accepted that radial pul- 
sations are the cause of the observed variations in 
many types of intrinsically variable stars. Model 
atmospheres have been calculated for stars earlier 
than the M spectral type. Dynamic computations 
require these static models as a basis from which 
to proceed. Radial pulsation theory has attempted, 
and succeeded quite admirably, to explain the 
observed variability of the hotter stars, especially 
the Cepheids. The literature is quite extensive on 
subjects such as period, stability of the pulsation, 
cause, amplitude and phase dependence of radius, 
temperature and luminosity. 

Ritter (1879) was the first to suggest that 
radial oscillations, with their associated surface 
temperature alterations, could provide the explanation 
necessary to describe the more regular periodic 
variables. Since that time, with a few roadblocks 
and detours along the way, and primarily due to 
the early work of Eddington (1917, 1918 a,b) and 
the later work of Zhevakin (1959, 1960) and Christy 
(1964) , the pulsation hypothesis has steadily become 
acceptable. Rosseland's (1949) book and recent 
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review articles by Zhevakin (1963) and Christy (1966) 
have summarized the pulsation theory of stellar 
variability. 

The earliest pulsational models assumed that 
the oscillation amplitude occurred throughout the 
entire stellar volume. The period was then pro- 
portional to p or f>^/2 m 1/2 * Fernie (1964) 

demonstrated that the period could be describee) fry 

. . 2 - 1/2 
being proportional to R M ' . This latter 

functional dependence could be explained by having 
the stellar atmosphere, and not the entire volume, 
participate in the oscillations. For the giant 
stars, the stellar core is inert. The pulsation 
amplitudes are damped and approach zero at the core- 
atmosphere boundary. In the Cepheids and Miras, 
although the partially degenerate core makes the 
extended atmosphere possible, the physical processes 
contributing to the variability occur in the stellar 
atmosphere . 

Many of the early models made approximations 
which, although allowing computational results, 
were not totally justified. It is now believe4 
that the radial oscillations of the stellar 
atmosphere are non-adiabatic and non-linear and 
may be solved only by modern computer methods . 

Christy (1964, 1968, 1969) attempted to soJ-VR 
the Cepheid variability problem with advanced 
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computers and techniques. His work, which was based 
on the contributions of many previous authors, appears 
to have been successful in explaining the cause, 
maintenance, and description of the pulsations. 

Christy's model treated non-adiabatic , non-linear 
motions of the stellar atmosphere which surrounds 
the stationary stellar core that supplies the energy 
generation. The observed temperature, radial 
velocity, luminosity and phase lag problems of 
Cepheid variables are adequately predicted by 
his results. 

As previously stated, theoretical model atmo- 
spheres for red giant stars have only recently been 
constructed. Without these static models for a 
basis, dynamical calculations have been almost 
non-existent. One of the first LPV models to 
include the ionization zones of HI and Hell and the 
dissociation of H 2 was constructed by Kamijo (1962) . 

His attempt to fabricate a red giant stellar model 
to explain Mira, itself, used linear, adiabatic 
pulsations and met with some success. 

Most recently, Keeley (1970a, 1970b) tried to 
improve upon the static and dynamic models by using 
non-linear, non-adiabatic equations, which were 
solved with the aid of a computer. He also made an 
approximate attempt to explain o Cet and was reasonably 
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successful. His model for Mira assumed a mass of 

one solar mass, an inert core comprising 54% of the 

3 7 

total mass, a luminosity of 10 ergs/second, a mean 
effective temperature of 2600 °K and a 300 day period 

To summarize Kelley's results: The radial 
velocity of the outer mass zone shells is character- 
ized by an asymmetric curve which has a rise time 
from maximum inward to maximum outward velocity of 
about 40% of a. cycle; the maximum rate or expansion 
occurs about 0.2 cycle after minimum radius; the 
maximum luminosity happens about 0.3 cycle after 
maximum compression, while minimum output occurs 
about 0.1 cycle before minimum radius; the luminosity 
curve has a steeper rise to maximum than decline 
to minimum and its range in bolometric magnitude is 
about 1™6. Unfortunately, Keeley was able to 
conclude little about the visual light curves or 
about the effective temperature and its changes. 

He stated that the phase relation of the luminosity 
to the radius is similar to that found for Cepheids. 

In order to compare quantities , table XXIV 
summarizes the results of Pettit and Nicholson 
(1933), the results of this study, and the predictions 
of Keeley (1970b) . 
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Table XXIV 

Intercomparison of Observational and Theoretical Results 

Pettit and 


Quantity 

Nicholson 

This Study 

Keeley 

Phase lag of 
maximum luminosity 
to minimum radius 

0.15 

0.22+0.03 

0.3 

Phase lag of 
minimum radius 
to minimum luminosity 

0.42 

0.32+0.04 

0.1 

AR/R 

0.37 

0.25+0.03 

0.45 

4m bol 

0™9 

0™7+0. 01 

1™6 

As evidenced by 

table XXIV, 

the agreement is 

quite good. Ho error 

estimations 

are listed 

in the 


table. Keeley's results were read from his figures. 

The radial velocity curves computed by Keeley 
appear similar to those of figure 15, although the 
agreement is far from excellent. In his models, 
shock waves were produced in the outer zones of the 
atmosphere. Shock waves have been observed in LPV 
stars (Deutsch and Merrill, 1959; Odgers and 
Kushwaha, 1960) and there is reason to believe that 
these shocks contribute to, or are the cause of, the 
high mass loss rates of about 2 x 10 Sl-./year, as 
calculated by Gehrz and Woolf (1971) . 



Both Keeley (1970b) and Christy (1968) give 
excellent descriptions of the physical processes 
which maintain the pulsations and produce the 
observed phase lags. It must be emphasized that, 
although the agreement is certainly not perfect, 
there appears to exist no major difficulty between 
observation and theory. The processes which result 
in variability seem to be similar for both the 
Cepheids and the long period variables. 



181 


XIV 

SUMMARY 

Infrared observations at an effective wavelength 
of 3 . 5y of a selected group of Long Period Variable 
(LPV) stars have been amassed throughout the full 
cycle of variability. Contrary to some beliefs of 
ten years ago, stars do vary in the infrared . Some 
of the LPV’s of this program altered their 3 . 5y flux 
by a factor of three. 

Although the radiation from Mira type variables 
did not repeat exactly in amplitude nor time of 
occurrence, the stars were sufficiently regular such 
that mean infrared light curves were produced. These 
mean curves were inter- compared for stars of spectral 
classification, M, S, and C. 

The 1 . 04y narrow band photometry of Lockwood 
(1972) and Lockwood and Wing (1971) was utilized to 
form a color index [1 . 04y ] - [3. 5y ] as a function of 
phase for stars which the two programs held in common. 
A color temperature scale was deduced from this color 
index on the assumption that a Lyra radiates as a 
10,000 °K black body at these two respective wave- 
lengths. The color temperature derived in this manner 
also assumed that the opacities of the two wavelengths 
were the same and that the cool star radiated as a 


black body. 
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From this color temperature through the Planck 
formula and from the 3.5vi magnitudes converted to 
absolute fluxes, the angular size of the source 
was calculated. With these two parameters, apparent 
visual and bolometric magnitudes as functions of 
phase were then calculated for the Mira variables. 

These quantities, as well as the computed radial 
velocities, were compared to observational quantities 
whenever possible. A color temperature-spectral 
type relation was also determined for the LPV stars 
and agreed adequately with the same relation as 
deduced for the non-variable or mildly variable M 
giants. The visual blanketing attributed to titanium 
oxide was also calculated as a function of spectral 
type and color temperature. 

The results for the Miras were as follows: The 
color temperature reached maximum simultaneously with 
maximum visual light {defined to be zero phase) , while 
minimum temperature was attained almost exactly one- 
half cycle later; the maximum compression (minimum 
radius) occurred about one-tenth cycle before maximum 
temperature, whereas maximum radius was about four- 
tenths of a cycle after minimum radius; maximum 
luminosity happened about one-tenth cycle after maximum 
temperature while the minimum energy output took place 
nearly one-half cycle later. The fractional change of 
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temperature was typically 20% while that of the 
diameter was somewhat greater and the luminosity 
altered by about a factor of two. 

These results are in excellent agreement with 
the thermocouple results of Pettit and Nicholson 
(1933) . The phase relation of maximum luminosity 
to minimum radius is similar to that observed for 
the extensively studied Cepheid type variables and 
possibly suggests that the physical mechanism pro- 
ducing the alterations might be similar. The cause 
of the variability, which must act like a Diesel 
engine absorbing heat upon compression and releasing 
heat upon expansion, appears to be Eddington's 
opacity valve. According to the models of Keeley 
(1970a, b) , this valve is the region of H and neutral 
hydrogen . 

Due to the general scarcity of models concerning 
long period variables , unfortunately little may be 
stated to support or denounce particular assumptions. 
It is hoped that this project will inspire more 
completely developed theoretical models which will 
explain the variations of these bloated stars. 
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XV 

SUGGESTIONS FOR FUTURE WORK 

By and large, the results presented here 
confirm those of Pettit and Nicholson's (1933) work on 
diameters and luminosities of long period variable 
stars. This present work was able to modify only 
a few of their results. 

The standard plea is for more observational data 
of all types throughout the entire cycle of variation. 
An improved method of temperature determination would, 
of course, be of the most practical use. A temperature 
scale must be developed for C and S stars so that 
physical parameters and variations may be compared 
to similar quantities for M stars. Spectra capable 
of radial velocity measurements are needed near 
minimum light where the available data are most scarce. 

There exists another experiment which could serve 
as a direct check on the validity of the assumptions 
which resulted in the computed angular diameters. A 
fifty foot interferometer should be able to resolve 
many of the objects listed. o Cet, R Cas, R Leo and 
X Cyg have diameters exceeding 25 milliseconds of arc. 
The interferometer could follow these variables 
throughout their respective cycles and could hopefully 



confirm at least the shape of , if not the absolute 
calibration of, the angular size versus phase 


curves . 
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